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1 Abstract 

Q ' It is shown that the Hamiltonian of the Einstein affine-metric (first order) formulation of General 

!>■ ■ 

Relativity (GR) leads to a constraint structure that allows the restoration of its unique gauge 
invariance, four-diffeomorphism, without the need of any field dependent redefinition of gauge 
parameters as is the case for the second order formulation. In the second order formulation of ADM 



5-H . 

gravity the need for such a redefinition is the result of the non-canonical change of variables [arXiv: 
0809.0097]. For the first order formulation, the necessity of such a redefinition "to correspond to 

> . 

' diffeomorphism invariance" (reported by Ghalati [arXiv: 0901.3344]) is just an artifact of using 

ON ' 

the Henneaux-Teitelboim-Zanelli ansatz [Nucl. Phys. B 332 (1990) 169], which is sensitive to the 
{Sj . choice of linear combination of tertiary constraints. This ansatz cannot be used as an algorithm 

Q> I for finding a gauge invariance, which is a unique property of a physical system, and it should 

o 

not be affected by different choices of linear combinations of non-primary first class constraints. 



C3 



^ ■ The algorithm of Castellani [Ann. Phys. 143 (1982) 357] is free from such a deficiency and it 

leads directly to four-diffeomorphism invariance for first, as well as for second order Hamiltonian 
formulations of GR. The distinct role of primary first class constraints, the effect of considering 
different linear combinations of constraints, the canonical transformations of phase-space variables, 
and their interplay are discussed in some detail for Hamiltonians of the second and first order 
formulations of metric GR. The first order formulation of Einstein-Cartan theory, which is the 
classical background of Loop Quantum Gravity, is also discussed. 
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I. INTRODUCTION 



We reconsider the Hamiltonian of General Relativity (GR) by using its equivalent first 
order form, the afhne- metric formulation of Einstein The reason for returning to this 
old and apparently solved problem is twofold. 

Firstly, in the literature this problem is claimed to have been solved 50 years ago by 
Arnowitt, Deser and Misner (ADM) [2J. The comparison of the ADM Hamiltonian with 
the Dirac Hamiltonian [3] shows some similarities (not equivalence) \&\; but Dirac based his 
derivation on second order, metric, GR. Preliminary results on the GR Hamiltonian for an 
equivalent first order formulation, based on slightly different but equivalent set of variables 



fa linear combination of affine connections 



the necessity to have tertiary constraints 



leads to a different conclusion; in particular, 
6|, contrary to the ADM treatment of the same 



problem. Recently, and for the first time, the Dirac analysis of the first order formulation 
was completed by Ghalati and McKeon |7| with an explicit demonstration of the closure of 
the Dirac procedure j8( and with the explicit form of the tertiary constraints given. This 
differs from the ADM Hamiltonian formulation and the reason for this discrepancy lies in 
the solving of the first class constraints, as indicated and discussed in {5, 6, 7|. According to 
the Dirac procedure, only second class constraints can be solved and the Poisson brackets 

nn 

of the remaining phase-space variables might be modified (Dirac brackets) [9|, [10J ■ 

Secondly, the Dirac analysis for systems with first class constraints cannot be consid- 
ered complete without the restoration of gauge transformations that, in accordance with the 

Dirac conjecture [8J, needs all first class constraints. Some steps of such a restoration and 

.00 



gg , was recently reported 



partial transformation for one phase-space variable, h 00 = 

by Ghalati 11] based on earlier obtained in [tJ first class constraints. In ll| the transfor- 
mation of h 00 was found using the approach proposed by Henneaux, Teitelboim and Zanelli 
(HTZ) 12J. This transformation is different from four-diffeomorphism and the author of 11] 
concluded that a field dependent redefinition of gauge parameters is needed for the derived 
transformation "to correspond to diffeomorphism invariance" This result is puzzling 

because in the Hamiltonian formulation of the second order GR, the necessity for such a 
field dependent redefinition is the result of a non-canonical change of variables |4|. Without 
such changes, the four-diffeomorphism follows without any redefinition of gauge parameters 
as was demonstrated for the Dirac formulation Q], as well as for the oldest Hamiltonian 
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formulation of GR due to Pirani, Schild and Skinner (PSS) |k| in [lllijj]. The equivalence 
of two formulations was demonstrated in [jjj] . In formulation of pi , the redefinitions of 
phase-space variables that might affect the result obtained in are canonical, so one ex- 
pects the complete and direct restoration of four-diffeomorphism without the need for any 
field dependent redefinition of gauge parameters, as is found in the second order formula- 
tion. The gauge invariance is a unique characteristic of a theory and equivalent second and 
first order formulations should give the same gauge invariance. The only difference in the 
second and first order formulations (considered in 0, [l^] and 11], respectively) lies in the 
methods of the restoration of gauge invariance that were used. For the Hamiltonian of the 
second order GR, the Castellani algorithm 16] was used for both Dirac 4] and PSS [bdT ] 



formulations; whereas in 



ll| the HTZ ansatz was employed. 



We want to clarify these two discrepancies simultaneously; and this dictates our choice 
of first order formulation. We perform our analysis for the standard and more familiar 
afline-metric formulation due to Einstein [l| because this formulation was the starting, point 
in the ADM analysis The affme-metric formulation, and one used in j^, [f], Q, 11]. 
are both equivalent to the second order formulation. The equivalence of afline-metric and 
metric GR was demonstrated by Einstein pj and for a different variable which is a linear 
combination of affine connections used in |7], the equivalence was explicitly demonstrated in 
Appendix A of [5|. Both first order formulations lead to similar Hamiltonians, as will become 
clear in the course of the calculations for the affme-metric formulation presented here and 
from the comparison of our results with ones obtained in [3]. Such a choice of variables 
(affine connections or a linear combination of affine connections) cannot be responsible for 
appearance of a different gauge invariance. There are some purely technical advantages in 
the parts of the calculations for one formulation over another that we will comment on in the 
course of our calculations; but they are not crucial, and neither formulation gives any overall 
advantage in calculation efficiency. The only "advantages" of affme-metric formulation that 
we want to mention, is the manifest covariant structure of primary first class constraints 
and possibility of performing a direct comparison with known transformations for affine 
connections. 

We perform this analysis for all dimensions higher than two as a specialization to the four 
dimensional case does not have any advantages or peculiarities. The standard expression 
"four-diffeomorphism" that we will use is equally well applied to all dimensions. The main 
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goal of our article is a thorough analysis of all the steps of calculation and, in particular, of the 
canonicity of all changes of phase-space variables that are used. All of these "technicalities" 
are our main concern. Our article is not an essay on the Hamiltonian formulation, but it 
is the Hamiltonian formulation itself, with all the steps of calculations and with sufficient 
details that anyone can repeat or check our derivations. So this article is the detailed 
proof that the Hamiltonian of the affine-metric GR and its constraints lead directly to four- 
diffeomorphism invariance without any redefinition of gauge parameters. This is exactly as 
it was shown in for the second order formulation of GR, if we use exactly the same 



algorithm of restoration of gauge invariance 



3- 



The plan of the paper is as follows. In the next Section, by using Dirac procedure, the 
Hamiltonian of the affine-metric formulation is obtained by performing the Hamiltonian 
reduction, i.e. the elimination of the phase-space variables associated with second class con- 
straints. In Section 3 closure of the Dirac procedure is demonstrated for a particular choice 
of tertiary constraints and the algebra of Poisson brackets (PBs) among tertiary constraints 
is compared with a similar algebra of secondary constraints of the second order formulation. 
In Section 4 the effect of different choices of linear combinations of non-primary 1 first class 
constraints are considered and their interplay with canonical transformations of phase-space 
variables is discussed. Section 5 provides arguments for a special role of primary first class 
constraints in the Hamiltonian formulation of gauge invariant theories. The examples from 
the Hamiltonian formulation of the second order metric GR as well as the Hamiltonian for- 
mulation of the first order Einstein-Cartan theory are discussed and the effect of unjustified 
manipulations with primary constraints is illustrated by an example of the Hamiltonian of 
oop Quantum Gravity (LQG). In Section 6 using the Castellani algorithm we restore as in 



ll| the partial transformations, but for all the phase-space variables of the reduced Hamil- 
tonian and we obtain all terms with temporal derivatives of fields and gauge parameters in 
the gauge transformations that coincide with four-diffeomorphism. We show that different 
choices of combinations of non-primary constraints do not affect the gauge transformations 
if we use the correct method to restore gauge invariance. In Section 7 we demonstrate the 
sensitivity of HTZ ansatz to a choice of tertiary constraints, contrary to the Castellani al- 



1 We call non-primary constraints all secondary, tertiary constraints, etc. In the literature they are some- 
times all called secondary; but because we consider two particular formulations with secondary and tertiary 
(not some general case) this terminology seems to be preferable to avoid confusion. 
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gorithm. Such an ambiguity and dependence of gauge invariance on a choice of a linear 
combination of non-primary first class constraints is in contradiction to the uniqueness of 
gauge invariance, which is an important property of a theory. This explains the origin of 



the puzzling result reported in The reason for the failure of the HTZ ansatz is dis- 

cussed, which is also related to a special role of the primary constraints - the true Masters 
of Hamiltonians for gauge invariant theories. In the Appendix A the details of the solution 
for the secondary second class constraints that were used in Section 2 to find the reduced 
Hamiltonian is given. 



II. THE HAMILTONIAN 



We start from the first order, affine-metric, Einstein action 



S E (g a ",r x aa )= jL(g^,T x aa )dx D (1) 
with the Lagrange density function 

l ( g ^, rt) = -V^gg aP (r^ A - + r* A r^ - t^t%) (2) 

where the metric, g^, and affine connection, T^, are treated as independent variables, 
g = det (g/ja/), and D is the dimension of spacetime. Greek letters are used for "spacetime" 
indices (/i = 0, 1, D — 1) and Latin letters for "space" indices {k — 1, D — 1). 

The first step in passing to the Hamiltonian formulation is the explicit separation of terms 
with temporal derivatives ("kinetic" part of L). For (T5]) we obtain 

J^kin — —V — 99 1 km,0 ~ V ~99 y- fc0,0 — 1 kmfi) + V ~99 1 Ofc,0- \ 6 ) 

This suggests the following field redefinition: 

= ^fcm, r° = 2Efc + r^, = Eg m — — — - 5^s o, (4) 

where £g m is a traceless field, = 0, i.e. E 00 = — r^. 

This redefinition does not affect the following components: 

-pfJ. _ -pM rm _ rm /e\ 
1 00 _ 1 00) 1 kp — 1 kp- \°) 
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After integration by parts and a change of variables, (jl]), the "kinetic" part of the La- 
grangian becomes diagonal 



L kin = {V^g km ) s fcm + 2 (V=gg k0 ) (0 s fc0 + {V=g~g 00 ) , ^oo = (v 73 ^) , *V (6) 

and "potential" part of the Lagrangian (terms without "velocities") is 



L pot = ( v^ 00 ) , fe r*, + 2 (v=*n , £ fc P - (v=^°) fc (r° - §-|e 00 ) 



(7) 



+ (V^~99 pq ) , rj, - 2 (v=^*) . (s 0p + r; m ) 



-gg km [TQ T, km + 2£ 0p r£ m + 2r^r| m - 4E 0fc E 



Orn 



D — 3 

/IV F 1 ? FP F9 OV V p _ _V V F p F<? 

Using the above variables and by performing Legendre transformation, we obtain the 
total Hamiltonian 



h t = g aP P aP + tapix* + r>nj° + t k 0m ul m + r%p% - (V^gg^) fi e«* - L pot . (8) 

The new set of independent variables, (HI) and (jSJ), and their conjugate momenta obey 
the fundamental Poisson brackets (PBs) 



{tf* (1?) , P, v (V)} = A$6 (1? - t) , IP"} = = - + , (9) 



{rgo, K } = {sL, n?} = *;a* - — l^&Ajg, {r£, n^} = s™a%. (io) 
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Here, only the first PB is written in a complete form. Further we will omit the delta 
functions and the dependence on "space" vectors x , y to shorten our notation, except in 
the cases where the derivatives of the delta functions appear; and where it is important to 
indicate with respect to what argument the differentiation is performed. 

We could equally well start from ([2]) to obtain the Hamiltonian and then perform canon- 
ical transformations in phase space from the original pairs {g a P,P a p) and (T^P" ") to a 
new set of variables: (g a ^P a p), (S Q/3 ,n^), (Sg m ,n° m ), (r£ ,nj°) and (r™,IT^). This 
is a canonical transformation that is automatically guaranteed for a linear and invertible 
redefinition of fields. 

As in any first order formulation, the Hamiltonian analysis leads to primary constraints 
equal in number to the number of independent fields 

Pa* - V^sjje^ « o, n Q ^o, nJ°^o, n° m ^o, n£«o (11) 

where we used 



— pa/3 uu 



(12) 



with 



E a/3 = A"' 3 



1 



af3 



(13) 



Among the primary constraints (fTTj) we have one pair of second class constraints with 
the following PB 



{Pa P - yf=gBQ^ n"} = -V=gE% , (14) 

which can be easily eliminated (they are of a special form; and therefore, the Dirac brackets 
among the remaining variables are the same as the corresponding PBs [10]). The solution 
for this pair is 



where 



W a = 0, 



'-9 



T&P p 

1 pLV 1 ' 



(15) 



j-a/3 — a a/3 
1 (MV — L -^\w 



D-2 



9^9 



a/3 ra/3 rpycr A7f 

1 \w — L ^ia>- 



(16) 
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Substitution of solution (fl5j) into the Hamiltonian ([8]) leads to (the first Hamiltonian 
reduction) 



Ht - ro n°° + Sg m n° m + r^n^ - L pot ^£ M „ - — =i^p a ^j . (17) 

The appearance of the combinations \/—gg alS and -j=I£P in the Hamiltonian and the 
invertability of I^f suggest the following canonical change of variables: 



h afi = V=ggr ,p i 7r, u = —I a JP aP . (18) 
Such a change obviously preserves the relation (using ffl2l and ffT6]) ) 

= <j^P Q/3 . (19) 



So, this transformation is canonical 13] and it can be also checked explicitly that 



{U*, 7^} = { v^^, -7= J ^ P p4 = A ^ ( 2 °) 



Equation ( 1191) is a necessary and sufficient condition for canonicity [171 ] . With such a 
canonical change, the Hamiltonian becomes much simpler; moreover, because the Lagrangian 
(T5]) is linear in y/—gg a ^, the Hamiltonian written in terms of new fields is polynomial. (So, 
;he often stated polynomiality of the tetrad Hamiltonian constraints in Ashtekar variables 
rs| as one of the advantages of the Hamiltonian formulation of the first order, tetrad-spin 
connection, Einstein-Cartan theory, is actually not something special since the first order 
affine-metric formulation is also polynomial in fields after the canonical transformation (jTSJ) 
is performed). 

Substitution of (fT8|) into (|Hj) constitutes the first reduction and gives the total Hamiltonian 
with fewer variables 

h t = r^ n°° + So m n° m + r^n^f + h c , (21) 

where the canonical Hamiltonian, H c , is 

H c = -r* (h™ + 2h 00 n k0 + 2h m0 n km ) - T° 00 (-/»<* + h 00 n 00 - h km 7i km ) (22) 
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-/^-^yvrooTToo " (>$ + h km TT km ) n 00 - (ih k \ 00 + Ah km n 0m - 2h%) n ok 



LOOymyA; , Q/,fcO /-pp ym pp ■spm\ riikrn snP OAP V 
Ok Om \ L mp^Ofc 1 fcm^OpJ — Z/1 ^kp^Qm ~ Zn ,k n 



Op 



+/i fcm (2r^r^ m — r^ p r^ 9 — T p k g r^ p ) — 2/i A:0 7r 00 r]™ n +2/i fcm (vr 0p r^ m — 27r ofc r^J —h p k T k pq +2h pk T 

After the elimination of the primary second class constraints (|T5l) (the first Hamiltonian 
reduction) and after performing the canonical transformation (fT8j) . we must continue the 
Dirac procedure and consider the time development of the remaining primary constraints. 
Two of them give the following secondary constraints: 

fC = {<, H c ) = -h ok + ^vroo - h km ir km ee (23) 

n° fe ° = {n° fc °, H c } = h™ + 2h°% 0k + 2h m % km = xl°; (24) 
and the only non-zero PB among them is 

{X° °,XT} = -XT. (25) 

The secondary constraints obviously (just from their field content) have zero PBs 
with all of the primary constraints: 

K°, K } = o, {xf, IT} = o, {xl°, n£} = o, (26) 

so they are first class, at least, at this stage of the Dirac procedure. 

Two remaining primary constraints, n° m ~ and IT^ « 0, lead to the secondary con- 
straints 

nr={nr,^e} = -^=xr, (27) 

nS = {nft.^c} = = xS- (28) 



pm' 



And, because H c has quadratic contributions in the corresponding coordinates (Eg m ,rj™), 
we have secondary constraints, which are second class, and two additional pairs of con- 
straints, (n° m ,x° m ) and ( n fc^Xm), which can be solved, and the corresponding pairs of 
canonical variables, (n° m , £q„J and (ilj™ , T™} , can be eliminated without affecting the 
PBs of remaining fields, as pairs (n° m , xl™) and v^kpi Xm) are a ^ so °f a special form. So- 
lution of these constraints is given in Appendix A. Substitution of the solutions into the 
Hamiltonian (the second Hamiltonian reduction) gives us the following total Hamiltonian 

Ht (C IT™ tf* vr Q/3 ) = + H c (29) 

with the canonical part H c 

H c = -TUT + (30) 
where H' c is the part of the canonical Hamiltonian after extracting the term — TqqX ® (sec- 
ondary constraints): 

Here we use the notation introduced by Dirac in the second order formulation of GR [3i] 
as it naturally arises in the course of calculation (see Appendix A): 

em = hn9 --fSxr> e m h qp = 6;. (32) 
Note that e nq is a short-hand notation, not a new variable. We introduce the combinations 
D® b and _D^ fe , and their explicit forms are calculated in Appendix and given by: 

D ok = hZ + h hk ^ (33) 



D» = -2h k H 0m + h% - j^D°J: - —D% + ^—j ( —D° c c - ttoo J (h 0q 5 k m + h k0 5l) . (34) 
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The next step of the Dirac procedure is the time development of the secondary first 
class constraints = {x°°;-^t}, to check whether they produce tertiary constraints. At 
this step of the analysis, we can proceed in two different ways. The direct way is to obtain 
-Hr} j and single out already known secondary constraints and form tertiary constraints 
from what is left. Another way is to isolate terms with secondary constraints, in addition to 
TqqX^ , in the total Hamiltonian before calculating Ht\- The terms in Ht proportional 
to the secondary constraints A^x°y will gi ye PB ^ooX^ } with the result which is 

proportional to the secondary constraints because of (|25|) . In this case, in order to find the 
tertiary constraints (if any), we have to consider the PBs of secondary constraints with the 
rest of the Hamiltonian, what is left after all terms with secondary constraints were isolated. 
Note that considering linear combinations of constraints is perfectly consistent with the 
Dirac procedure (as a linear combination of constraints is also a constraint). 

Performing the direct calculation, we use the following simple PBs of the secondary 
constraint Xo° with the combinations presented in (l30l)-(!34l): 

{*" £5o} = M ,* 00 } = MV°o} = -oo, (35) 

{x$ > ,h ah } = 0, {tf,I%} = 0, (36) 

{X™ Ah ko ir 00 + 4h km n 0rn - 2h pk } = Ah ko n 00 + Ah km n 0m - 2h pk , (37) 

as well as 

{X° °, h% - 2h k H 0m } = h% - 2h k ^ 0m (38) 
that, in combination with ( 1351) and ( 1361) . gives 

{X%°,D%} = D%. (39) 

Using definition (152]) we also find: 

{Xo°, e km } = e fcm , {x° °, h km ) = -h km . (40) 
The first PB in ( HOI requires calculations; but the second is just the result of e km h kn = 5™. 
Similar PBs for x k ° are: 

11 



{X° b °,e km } = 0, {xl°,h km }=0, (41) 



{xl°,h™}=0. (42) 
We would like to mention here that we originally analyzed the first order of the EH action 

fin 

in different variables [5|, |6( where we used another combination 2 

& = - \ {6$%, + tfrj) (43) 

and discussed the unavoidable appearance of tertiary constraints. These variables, (1431) . 
simplify the calculations because the "kinetic" part of the Lagrangian becomes diagonal 
automatically: Lkm = —\ / ~99 al3 £af3,o- But this simplification appears only in the first steps 
of the Hamiltonian analysis. These variables were used by Ghalati and McKeon to find the 
Hamiltonian for the first order formulation of GR and, for the first time, they explicitly 
demonstrated the closure of the Dirac procedure [3]. Firstly, our choice to start from the 
afline connections that lead to combinations (j3J) and (jSJ) is dictated by our goal to compare 
our results with [2j , where the authors claimed that they used Palatini formulation, 3 which is 
the affine-metric formulation due to Einstein. Secondly, the variables (HI)-© also diagonalize 
the "kinetic" part of the Lagrangian (J6j); and, what is more important, the primary first class 
constraints n°° appear in the covariant form. In addition, if we used £*/3 from (j4"3"l) instead 
of r^g then some of the brackets fl35l) - fl39l) would be more complicated. For example, fl39|) 
became non-local (proportional to derivatives of delta functions). In this case, it would be 
impossible to use simple associative properties of PBs for the terms of H' c in the second and 
third lines of (l3~il when calculating Whereas, using the Hamiltonian in terms of 

T^p and properties of the above PBs the calculation of {xo°, H' c } is greatly simplified, for 
example, from ([H and flUJ it immediately follows that {xo°, Dfh hz Df) = Dfh hz Df , etc. 
One more advantage in using the original variables is the simplification of restoration of 
gauge invariance; we will not need to restore the gauge transformations of from that of 



2 This is generalization of change of variables that we found considering the Hamiltonian of 2D in ■ Later 
we learned that such a combination was known before (see Horava 20( for any dimension and Kijowski 
for four-dimesional case 0). 

3 This formulation was originally introduced by Einstein but continues to be mistakenly attributed to 



Palatini [22| (see also Palatini's original paper [23| and its English translation [24[) 
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£a8- The above arguments are mainly technical and the real advantage is just a manifestly 
covariant form of the primary first class constraints. As we worked with both and T^g , 
we have to admit that there is no overall advantage in a particular choice of variables and 
the total amount of calculation (difficulties) is conserved. 
For the time development of Xo° we have 

X°o° = M°, H T } = -TUl°} + {xg°, K) = r k 00 xT + {x°o°, K) (44) 
where using properties fl35|) - P0|) . the last PB can be just read off 

{X° °, K) =H' C + {2h k \ 00 + 2h km n 0m - h^) k . (45) 

Note that the right-hand side of ( )45l) is not the Hamiltonian density and we cannot neglect 
this spatial derivative. 

Because all PBs ( |35|) -( 140|) are local (no derivatives of delta functions), using them and 
associative properties of PB makes the result f|45l) almost obvious. The first term in f|44l) . 
TooXfc ? is proportional to the secondary constraints. The rest is given by (1431) and to find 
out whether we have closure of Dirac procedure at this stage or if the next generation 
of constraints appears, we first have to find in (T4"5l) the combinations proportional to the 
secondary constraints. If there are contributions, which are not proportional to secondary 
constraints in (j4"5l) . then we have tertiary constraints. 

We proceed as follows. Using secondary constraints (|23l) and (|24l) we express some fields in 
terms of constraints and remaining variables, substitute them into the Hamiltonian, isolate 
the contributions, which are proportional to secondary first class constraints, and then work 
with what is left. Note that according to the Dirac procedure we cannot solve first class 



constraints, as it was done in [2| and which is shown especially clearly 



DV Faddeev in 



25|; 



These re- 



and the consequence of solving first class constraints was discussed in 
expressions, to single out contributions proportional to secondary first class constraints, after 
calculation of PB (f4"5l) or before, as it was done in Q, is consistent with the possibility of 
using linear combinations of non-primary first class constraints (the role primary first class 
constraints will be discussed in Section 5). 

The only two phase-space variables that can be unambiguously re-expressed (by algebraic 
operations) using secondary constraints fl23|) . (1241) are two momenta: 



13 



(46) 



^ ' ^00 1 imO^ , ■'■ 1 ,,00 



(47) 



After substitution of ()46l) and (I4T1) into the Hamiltonian (l3il . we separate terms pro- 
portional to an d write the rest of the Hamiltonian as a sum of three contributions of 
different order in it^m (e.g. H" (2) is of the second order in 7Tk m , etc.). Performing the same 
operations with (f45j) we obtain the following form of the Hamiltonian H' c 



H> = H" (2) + H" (1) + H" (0) + AX° + ST T 
and equation (T4"5l) can be rewritten as 



(48) 



{xo 10 , K) = K (2) + K (i) + K (o) + {x° \ a^} 



(49) 



where 



(50) 



2 (e m % m<? ) k - 2 (e m % fc9 ) 



(e np h 00 ) 



(51) 



^(0) = ^ 



(52) 



;h ay h bz e™ {e a V) n —{ey^) x + 



1 1 



nx h ab (e«*/i») VtSoC^ 00 ) 



4 L> - 2 



.4° 



D - 1 



D-2Df h k0 h™ h k0 h m0 
00 D - 1 h 00 h 00 h 00 h 00 h 00 



(53) 



D 
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1 1 



D — lh 00 



{h 0a 5 b n + h"%) 



I PS nz 

^ h bzD a 



1 w»* 1 1 

-h ay h bz e nx D x +4^3^' 



h a bh yz D n 



A ~ 2 /^oo 71-00 2 /^oo ^Q" 1 + /^oo + /^oo /^oo + 2 /^OO /^OO 



(54) 



1 

2^ 



1 



e nk h ab D™h yz 



ab i 



AD -2 



and 



— J 

4 



-^£> fl -~h ay h bz e kx D x + -—— e kx h ab h yz D x 



£) = h kq —h 00 + 2h kq —h p0 ir + h kq - — D ok - — D 0q 



A and A k are the functions which depend on the fields and derivatives and their explicit form 
is not needed for further calculations in this article; they are given only for completeness. Our 
main interest is in the contributions, which are not proportional to secondary constraints, 
H" (i). The explicit form of S m is found by comparison of the parts not proportional to the 
secondary constraints in (|45|) with H' c 



S 1 ' 



-2—* kq 



^oo e nk i ~t~ z f^o e nk i '",P fjOO 1 ">P ^00 1 ^00 V- ),k 
And, as in the previous steps of the Hamiltonian reduction, writing h km in terms of e km 
makes the expressions more transparent. Note that (148]) . contrary to (1451) . is the Hamiltonian 
density and a surface term can be neglected in subsequent calculations. 

Using (|48j) and (!50|) - (l52|) the calculation of the time development of x k ° is straightforward 



^ + i>7m + m(e km in, l (55) 



where 



oo 



(56) 



{xT,H' c }=2(e mq 7r. 



mqj k 



- 2 (e mq 7T kq ) 



h 00 



(e^ 00 ) fc7 r np +{ X ; ,^ }. (57) 



Note that {x° k ° , H' c (2)} = { X ° k °, H' c (0)} = which is based on simple PBs of xl° with 
combinations of fields presented in (j4T]) - (jl2|) . 
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The last term in (1571) gives contributions proportional to whereas the first three terms 
in fl57|) cannot be expressed as a linear combination of moreover they coincide with the 
expression in square brackets of H" (1) in ( 15TT) . This part of ( 1571) . which is not proportional 
to the secondary constraints, can be chosen to be called a tertiary constraint 

r°° 2 {e m H mq ) M - 2 {e m H kq ) m - ± (e"^°%7r np . (58) 
Taking into account (|58|) . the bracket {xo°,H' c } can rewritten as 

M°, K) = K (2) + ^rf + H>> (0) + {*» A^} , 

so the terms which are not proportional to already known constraints should be called a 
new, tertiary constraint. Using the exact expressions of H" (2) and H" (0) from (1501) and 
(152]) we can name the following combination as a tertiary constraint 

r oo _ e k Pe m q {vkq _ Vmq) _ (e km h 00 ) mk + \ (eV) y fc, JL (e^h 00 ) x (59) 



--h ay h bz e nx (e ab h 00 ) ^ (e^h 00 ) + --^— e nx h ab (e ab h 00 ) h yz -^ (e**h 00 ) . 
^ a v oz \ / ,n ^oo v < x A D — 2 ' n hP° > x 

The canonical Hamiltonian written in terms of the constraints and r°° is 

H c = -Cx™ + ^ + + A"*" (60) 

Note that the choice of tertiary constraints is not unique. For example, if we start 
from Xo° we would name the whole combination H" (2) + H" (1) + H" (0) in ( )49l) a tertiary 
constraint Tq°, because in this case r°° has not yet been found. Such arbitrariness in the 
choice of constraints looks ambiguous. Firstly, it does not contradict the Dirac procedure 
as any linear combination of constraints is also a constraint. Secondly, and we will show 
this below, if the correct method of the restoration of gauge invariance is used then the final 
result does not depend on a choice of tertiary constraints. 

We have chosen such a form of tertiary constraints, (15"8~1) and (1591) . because calculations 
of PBs among secondary and these tertiary constraints is almost manifest due to the simple 
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properties of the PBs of the secondary constraints with their combinations presented here 
(e.g. (j35])-(j40D and (p]) - (p]l ). It is easy to show that 

«V°°}=0. (61) 

The brackets {t°°, A^\ ^ } are proportional to the secondary constraints because of ( 16T1) . 
In addition, because of the relatively simple form of r°° ( fl58l) and (159]) ). the calculation of 
PBs among them is not inordinately tedious. We must find these PBs to prove the closure of 
the Dirac procedure; and if it closes for one choice of tertiary constraints, then it closes for 
any combination of them. For these constraints, r°°, the only possibility to have a closure 
is to demonstrate that 



KV°°} = or ~r°°, (62) 

as PBs of combinations (!58l) and (|59l) cannot form secondary constraints just because of 
their field content. In the next Section, we consider the calculations of PBs (162]) and the 
closure of the Dirac procedure. 



III. CALCULATION OF {t°°,t°°}> CLOSURE OF THE DIRAC PROCEDURE 
AND ALGEBRA OF TERTIARY CONSTRAINTS 



Even considering relatively simple combinations r° ( ( 1581) and ( 159]) ) as a choice of tertiary 
constraints, the calculation of the PBs among them is a laborious procedure. This fact and 
a variety of other choices of tertiary constraints were the reason why the first attempts to 



orove a closure of the Dirac procedure for first order formulation of GR were not 



, |6j where the variables 



finished 
2611 and 



were used; only later the proof was completed in j?]. 
the expectations outlined in 0, (| are thus realized. The presence of derivatives of fields 
imposes some additional complications and the best way of dealing with such calculations is 
to use test functions [9| that were demonstrated in some detail for constraints of Yang-Mills 



theory 27|. More details of calculation using test functions were given and applied to the 
first order formulation of GR in [3, [ll|. The PBs among constraints can be written in the 



form with the explicit presence of test functions (e.g. see work of Faddeev |25j, Section 3), 
which might be useful for further calculations (e.g. restoration of gauge invariance, analysis 
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of different choices of tertiary constraints, etc.) compared with the standard form which 
contains derivatives of delta functions. 

Introducing test functions, f(x) and g(y), that have a zero PB with phase space 
variables, we calculate 



<fi?dif {/ (1?) r" , g (V) r°° (V) } 





dxdy (...)-£_> 5 (a; - y) = j dx (...)^ (63) 
where, to shorten the notation, we omit the integrals, i.e. 

Oft = (■»)* . (64) 



After long but straightforward calculation, f[63|) can be presented in the following form 
for //, v = 0, z: 

{/ (!?) r°° (1?) , 5 (V) Tf (V)} = fjgT? - fifg 4 , (65) 
{/(l?)r 00 (-?),^(^)r 00 (^)} =/feV;-/, t WV; , (66) 
{/(l?)r 00 (l?)^(^)r° (^)} = - /^r 00 , (67) 

{/ (1?) r°° (1?) , 9 (V) r 00 (V) } = /,^r 00 - fg^° . (68) 

From these expressions we can also obtain the standard form of PBs with delta functions 
if we rearrange (1631) in the form without derivatives of the test functions (note that to do 
this for ([63|) - ([5S]1 derivatives of delta functions are unavoidable) 

J J ' <n?<nt{f{-$)Tf{-3) ,gCy)r° Ot)} = I Jd2djtf(l?)(...)- Stit g(jr). (69) 
For ([SSD-dSTD we obtain: 

K° ("?) , r°° (V)} = (x) d»5 (x-y)- r°° (y) d*5 (x - y) , (70) 
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{r 00 ("?) , r 00 (V)} = -/> 00 (x) e fcp (x) r° p ° (x) 0*6 (x - y) + h 00 (y) e kp (y) r° p ° (y) d x k 5 (x - y) , 

(71) 

N° ( * ) , r°° (V)} = r 00 (*) fl?* tx-y)- r 00 (y) 9*5 (a; - y) , (72) 

K° {^) , r 00 (V)} = r 00 (x) flftf (x - y) - r 00 (y) 9f5 (x - y) (73) 

(here we use the notation: d% = or df = ^pr). As a consistency check we can integrate 
the above expressions with J d~y f (1?) g (y) which leads us back to the previous 

form (!65l) -( j67j) . So, these, (J65l) -( j68l) and ( J70l) - (|73l are two different but equivalent forms of 
the constraint algebra. 

This algebra of constraints (in one form or another) is equivalent with the algebra found 
in which becomes clear after canonical transformations are performed (see next Section) 
and it is also the same as the algebra of constraints (but secondary) given by Faddeev 25] 
and Teitelboim [28] (despite different expressions for constraints themselves). 

The conventional form of the algebra of secondary constraints for the conventional Hamil- 
tonian formulation of the second order EH action, known also as "Dirac's algebra" or "hy- 
persurface deformation algebra", is: 

{H L (x) , H L (x')} = e rs (x) H s (x) 5 >r[x) (x - x') - e rs (x) H s (x') 5, r( , (x - x') , (74) 



{H s (x) , H L (x')} = H L (x) 6 Ax) (x - x') , (75) 



{H r (x) , H s (x')} = H s (x) 5 tr{x ) (x - x') - H r (x) 6 >s (x>) (x - x) ; (76) 

where Hl and H s are the "Hamiltonian" and "diffeomorphism" constraints, respectively. 



The a 



32 



33]. 



gebra (174j) - (!T6"I) can be found in slightly different forms, for example, in [8j, 



29 



30 



31 



Comparing fl74j) - fl76]) with (J70ll - (j73|) one can notice a difference. It is not related to the 
variety of notations used in the literature, but to the obviously less symmetric form of (IT5]) 
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in the conventional algebra, contrary to (!73|) and (I68j) . This discrepancy must be clarified. 
Let us trace out the origin of the algebra (1T41) - (1T6"1) starting from its first name "Dirac's 
algebra" . It appeared for the first time in the Dirac book |8j , where he referred to his paper 



34J . But in [34] he derived this algebra for the motion of space-like surfaces, as PBs among 



tangential and normal to a surface variables, not as PBs among the secondary constrains of 
GR. Kuchar in 30], by geometrical reasoning, showed how 7Yl and 7i s ( "super-Hamiltonian" 



and "sup er- moment a" , in his terminology) "represent the set of deformations of space- 
like hypersurfaces" . Probably, after that the name "hypersurface deformation algebra" 
appeared. In 32] Teitelboim reconstructed this algebra by "a simple geometrical argument 
based exclusively on the path independence of the dynamical evolution" , i.e. on "the 'motion' 
of a three-dimensional cut in a four- dimensional manifold of hyperbolic signature" . In his 
later work 28] this algebra was altered by another one where the PB (1751) was replaced 



by (I7T1) . This transition was not explained and left unnoticed, which is strange, especially 
because when presenting the new algebra he referred to his old paper 32] where the algebra 
is different. This algebra of PBs among secondary constraints was derived by Faddeev in 
251 ] where he considered the first order formulation of GR; in his paper, the algebra among 
constraints is written in the form of fl65|) - fl68|) . The same algebra was presented by Ghalati 
and McKeon [7|. Our calculation also results in (165]) - (168]) . or equivalently in (170|) - (l73|) . 

In 35] analyzing the Dirac derivation of 34J, the author made a conclusion that "Dirac's 
derivation of the constraint algebra cannot be considered satisfactory" . This statement must 
be clarified. Dirac did not derive the algebra of constraints as they were not even known and 

n 

appeared only a few years after [3] ; and in this article there is no statement that constraints 
satisfy "Dirac algebra" of 3^|. This conclusion was made by other authors and without 
calculation. 

To answer the question why the discrepancy, the difference between ( 173!) and (1751) . ap- 
pears, we are planning to revisit our analysis of the Dirac formulation of the second order 
of EH action given in Q], where a different choice of secondary constraints was used. The 
results will be reported elsewhere. 

If we are interested only in a demonstration of closure of the Dirac procedure and in 
restoration of four-diffeomorphism, we need to consider the time development of all first 
class constraints started from primary, i.e. to calculate PB with Hamiltonian. In this case, 
the algebra of particularly chosen tertiary constraints is not important. In addition, one 
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can completely avoid the non-locality of these PBs (i.e. derivatives of test functions in 
fl65|) - fl68|) or derivatives of delta functions in fl701) - fl73l) ) as the PB of a constraint with the 
Hamiltonian is always local (one integration has to be performed as Hamiltonian in field 
theories is integral of Hamiltonian density). For example, PBs between tertiary constraints 
and the Hamiltonian is defined as 



. oo 



W°.^} = {/ 



As the Hamiltonian for the EH action is a linear combination of the constraints, what we 
actually need to calculate are the following PBs 



/ (?) (1?) , J dy>9» (V) r°° (t) I (77) 

where, according to our choice of constraints (I60I) . we have to put g°(y) = ^7=^ anc ^ 
9 (y) — , 00 /l>\ ■ Using (170|) - (!73|) we can easily find: 

>(*), / - - (^) , (78) 



r„ M (1? 



<fiW (y) s4j } = * + h °° ek ' T "' (79) 



AY, 



rra). I .UrT(7)^—^[^) ( -[f^) (SO) 

'Hi r , / JT?~00 /"r+\ ^ 1 . _00 | M /^00 1 



rfjtf m _J=-^_j r ^_^ . (81) 

Equally well we can use the PBs, (|65l) - (l68l) . where by integration by parts we have to 
move the derivative from a test function /, and put / = 1 at the end of the calculations (see, 
e.g. 13 )• This will result in the same expressions fl78|) - flHT]) . We also want to emphasize 
that these calculations can be performed directly, without any reference to a particular form 
of the algebra of constraints, which depends on our choice of constraints (see next Section). 

The above PBs, (1751) - (15T1) . complete the proof of closure of the Dirac procedure that the 
PBs of the tertiary constraints with the Hamiltonian are proportional to already known 
constraints 
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K°. Br) = {r™, ^r»» } + {r» ^r»} + {r» A"} X „ . (82) 

In next two Sections we will discuss the role of different choices of constraints and canoni- 
cal transformations, which is important in general, but also will be needed for the discussion 
of methods of restoration of gauge symmetry (Sections 6 and 7) that, in accordance with 
the Dirac conjecture |8|, is generated by the full set of first class constraints. 

Both equivalent forms of the PB algebra among tertiary constraints, which we chose out of 
many possible combinations, (|65l) - (l68l) and ( J70l) - (173l) . might be useful in calculations that we 
need to perform, especially for the general analysis of a role of different linear combinations 
of non-primary first class constraints. A particular choice of tertiary constraints can lead to 
considerable simplification in some parts of the analysis; but the algebra with derivatives of 
delta functions or test functions can be completely avoided j^, [l4], and are not needed for 
oroof of closure nor for the restoration of gauge invariance using the Castellani algorithm 
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IV. LINEAR COMBINATIONS OF TERTIARY CONSTRAINTS, CANONICAL 
TRANSFORMATIONS AND THEIR INTERPLAY 

Considering the time development of secondary constraints in Section 2 we demonstrated 
that tertiary constraints can be defined in different ways and all such combinations do not 
contradict the Dirac procedure. All the different choices are linear combinations of each 
other, i.e. the Hamiltonian formulation of affine-metric GR provides an example of the 
theory with the non-artificial appearance of different linear combinations of constraints. We 
would like to discuss this apparent ambiguity of the Hamiltonian procedure. We restrict our 
discussion to different combinations of tertiary constraints that we initially defined in (1581) . 
( |59|) and the corresponding part of the canonical Hamiltonian 

1 h ok 

Till _ 1 ^-00 , " _00 /QQ\ 

H c - ^0 T + ^0 T k ■ ( 83 ) 

Of course, more choices exist if the secondary constraints are also used in such redefini- 
tions; but it will just make the calculations more involved. Some conclusions can be made 
based on simple examples. Note that the role of first class primary constraints is quite 
special and will be discussed in next Section. 
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One possible choice of tertiary constraints is 



r 00 - J_r 00 r 00 - r 00 (M) 
In terms of these constraints, the part, H", of the canonical Hamiltonian (1831) is 

LOfc 

^ = -o°° + W , Tf - (85) 

The PBs among constraints f and can be easily found by using their relations with 
the original choice and the corresponding algebra of constraints (!65|) - (l68|) . e.g. 

{/(l?)ro(-?),^(V)r° (V)} = -/^°Vo 00 + /,^ 00 - /f °V, • (86) 

The whole algebra of Tq° and f°° can be calculated and the closure of the Dirac procedure 
can be demonstrated using these combinations. 
Another choice is: 

1 h ok 

^00 _ 1 _00 , n _00 -00 _ _00 /o 7 N 

that leads to a very simple expression for the corresponding part of the canonical Hamiltonian 



H>> = f 00 , (88) 

and again the PB algebra of constraints Tq° and f°° can be found using (!H7|) and (I65l) -(l68l). 
e.g. 

{/(l?)f 00 (l?),^(V)rr(V)} = 

uOk uOk UOk LOO 

f ~00 -f^00„ I £ n ~00„ ^" ~00„ , f' ,i ~00 x„ >« ~00 / QCl \ 

J,i9T - /r + fj^n g ti - f—n g, k + I9j^r k - fg— TQ . (89) 

Is there any physical significance in a particular choice of constraints? Of course, there 
are some possible technical (computational) advantages; but considering different linear 
combinations of tertiary constraints should not affect the physical results. The simplest 
argument is to convert our reduced total Hamiltonian into the corresponding Lagrangian by 
inverse Legendre transformation 
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L = n Il™ + h^ap -H T = h^ap - r« oX °° - H" - A^x™, (90) 

where H" is a functional, H" (/i 00 , h ok , h km , 7Tfc m ). Whatever combination we consider, fl83l) . 
(!85|) or (188|) . we have the same Lagrangian (H" = H" = H") and by calling some parts of 
these Lagrangians by 'tertiary constraints', which is nothing more than a short notation at 
the Lagrangian level, we cannot affect the physics; and in particular, the gauge invariance 
should not change. So, if a choice of tertiary constraints cannot influence gauge invariance, 
then an algebra of PBs for this particular choice should not bear any physical significance. 
For the Hamiltonian formulation of the second order metric GR, we can make the same 
conclusion, but, of course, for possible choices of secondary first class constraints and their 
algebra, contrary to a broadly accepted view that a particular choice of constraints/algebra 
has some physical significance that is even reflected in special names given to one particular 
choice, "Hamiltonian" and "spatial diffeomorphism" constraints. These particular combina- 
tions become special only after a non-canonical change of variables is made (see discussion 
in next Section). Related to the combinations of constraints idea of the Master Constraint 
Programme is of limited interest as any physical results cannot depend on a particular choice 
of non-primary first class constraints; in other words, in the "society" of non-primary con- 
straints there is no place for a Master. One additional conclusion that is connected to the 
freedom to choose combinations of tertiary constraints is more technical and related to the 
methods of the restoration of gauge symmetry, based on a full set of first class constraints 
(the Dirac conjecture). Gauge invariance should be independent of a choice of non-primary 
first class constraints. So any method, which is sensitive to redefinition of constraints (gives 
different transformations) is not correct (see, Sections 6, 7). 

For any Hamiltonian formulation, if change of phase-space variables is performed then it 
must be canonical. Also such a change should not affect the physical properties of a system, 
in particular, it should not change its gauge invariance. In the second order formulation of 
GR we considered the connection of two Hamiltonians, Dirac's and PSS, and demonstrated 
that they are related by a canonical transformation and both lead to the same gauge trans- 
formation which is four diffeomorphism invariance [4, 14, 15|. In the passage from PSS to the 
Dirac formulation, we followed Dirac's idea: to simplify primary constraints. He achieved 
this by modifying the ori gina l Lagrangian. We worked in phase space and performed the 
canonical transformation 15|. In first order formulation of GR that we consider here, the 
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primary first class constraints already have the simplest possible form: pure momenta con- 
jugate to Tq , and any further simplification is impossible. We should have different reasons 
to look for canonical transformations. One such a reason is to simplify the expressions for 
the secondary constraints, for example: 

,00 LOfc , 7 OCI h km ^ hOk j i 00- -00 / m ^ 

Xo = ~h,k + h ^oo - n 7r fcm = -h >k + h 7r 00 = X > (91) 

,00 lOO | qlOO^ , numO^ 7 00 , r>7 00- -00 /r,0\ 

i.e. to introduce new momenta 

Lfcm 

7T00 = TTOO - -J^0 n km , (93) 

h m0 

TTo/c = T^Ok + -J^^km ■ (94) 

These two redefinitions, fl93l) and ([M]) . are algebraic and invertible; but this is not enough 
to preserve the canonicity of phase-space variables, and it must be accompanied by a change 
of the remaining phase-space variables, which are involved in redefinitions ( l93l -( |94l) . Such 
a necessary and sufficient condition for the transformation to be canonical is 17] 



ti$7r a i3 = h a p n al3 . (95) 
Let us restrict our search for canonical transformations by assuming that 

h 00 = h 00 , h ok = h ok . (96) 

(This restriction has to be relaxed if it is not possible to satisfy ( |95l) .) After substitution of 
new variables in terms of old into (195]) . and some simple rearrangements we obtain 

h k ^ km = (h 00 h km - h ok h 0m ) ^n km . (97) 
To fulfill this condition we have to define: 

jr km = -^vr fcm , h km = h 00 h km - h ok h 0m = h 00 e km . (98) 
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Of course, that the transformation which involves (1931 . f|94|) . (1961) and (|98|) is canonical 
can be checked by direct calculation of the PBs among all new variables; and this should 
lead to: 



h^,n u A = \h^(h^,7r^),^(h al3 ,7r u M ={^,7^}, (99) 



^} = j>/? ^00 ) 



= o, 



{tTq/3, TT^m) = i 71 "^ (^ Q/3 ' 7r ^) ' ^ 7r ^)} /l a0 ^ = 0. 



Note, that the second equation in (l98i) was used by Faddeev |25J], but without any discus- 
sion of canonicity, i.e. necessary changes for the rest of the variables. In Faddeev's approach 
all of the problems related to such a change of variables are hidden because, in addition, 
some first class constraints were solved. This is against the Dirac procedure and it eliminates 
the possibility to restore gauge invariance as all first class constraints are needed (a simple 
example can be found in [51]). 

Equally well, instead of a simplification of secondary constraints, we can try to use a 
combination e km which naturally appeared when the secondary second class constraints 
were solved (the same combination, e fcm , was used in the second order GR by Dirac [3]). We 
used e hm as a short-hand notation; but it is possible to find such a canonical transformation 
that converts e km into a new variable. Let us introduce a variable 

h km = e km = h km __ ^qq) 

As in the previous case, we will restrict our search by imposing 

h 00 = h 00 , h ok = h oh (101) 

and use the same condition as (1951) 

h°o7r a/ 3 = ■ (102) 

Substitution of h km from (HOOD into fTT02]) gives 
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h k ™TT km + 2h°Q I TT k + -J^T^km J + tifl ( ^00 ~ ^ 00 ^00 71 ' km ) = ^O™^ + 2 ^°0%fc + /i° °7foO 

(103) 

and the redefinition of momenta follows: 

TTfcm = TTfcm , (104) 
TTOfe = + J^^km , (105) 

= - • ( 106 ) 

It is easy to check that fll04p - (1106p . together with ( 11001) and f llOll) . give the same PBs as 
( 1991) and so these transformations are canonical. 

Another possible argument to find canonical transformations is to look at expressions 
for tertiary constraints (or rather our first choice of tertiary constraints) (1581) . (1591) and try 
to simplify them. There is one obvious combination, h 00 e km . And it can be used to build 
canonical transformations; but it will lead to the same transformation that we have already 
considered in the first example (198|) . If we use it together with (I93 |) -( l94|) . we will get the 
canonical transformations, which simplify secondary and tertiary constraints simultaneously. 

Substitution of (1981) into constraints (158]) and (1591) gives: 

f 00 = 2 f~ h m q ~ m \ _ 2 (^J -hjlt np , (107) 
V / ,k V / ,m 

T~q — H P H 9 ij^mp'^kq ^kp^rnq) 



- htm + \^kJ% - h ay hJi nx hih% + ^r^ftjl^ 2 (106 

where hkm is defined as h gm h mp = 5 g and related to the original variables by hkm — jjfttr-. 



These constraints, ( I107p and ( I108p . are very similar to the constraints obtained in 11] 
for the first order formulation based on variables (T4*3i) . This is not a surprise, as both first 
order formulations are equivalent to the EH action. Up to a simple rearrangement and with 
a different notation the expressions corresponding to fllQ7p -f lT08l) in [ll| are the following: 
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(11071) is equivalent with Eq. (56) of [11 



the corresponding term in Eq. (58) of 



and in (11081) . only the first term differs in sign from 
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The effect of canonical transformations on the algebra of constraints is simple and it 
preserves its form (form-invariance). In general, a canonical change of variables leads to 
changes in constraints 

(q, p) -> (Q, P) : ipfj, (q, p) —>■ ^ (Q, P) (109) 
and if the algebra of constraints in old variables is 

{VV> Vv} = c^i/jj , (110) 

where c 1 are structure functions, then in new variables its form should be preserved and 
given by 

= C£,¥ 7 (111) 
with the simple condition on the structure functions 

£>/ (Qi P) = c ]iv (<liP)q=q(Q,P), p =p(Q,P) ■ ( 112 ) 

For the first time such properties were demonstrated for the canonical transformation 



in linearized gravity 



36j and later for a complete formulation [15j. In the Hamiltonian 
formulation of the first order EH action, the form-invariance of the algebra of constraints 
is also preserved after canonical transformations. For example, our secondary (1231) - (12^1) 
and tertiary r°° (|58l) -( l59l constraints have the algebra of PBs given in (]2"5]1 . ( 16T1) . ( I7D1 - II731 . 
Performing the canonical transformations ( 1931) . ( |94|) . ( l96i) . ( 1981) . we obtain new secondary 
fl9Tj) - fl92l and tertiary f°° (I107p - (I108I) constraints and calculate new PBs among them. 
The algebra of new PBs is related to the old one exactly as is described by fl 1 09 [) - (1 1 1 21) . i.e. 
it is form-invariant. 

In the first order formulation of the EH action, we consider the canonical change of 
variables that does not affect the primary variables - fields for which their corresponding 
momenta are primary first class constraints. So the corresponding algebra of all first class 
primary constraints is form-invariant (in new and old variables primary first class constraints 
have zero PBs with the rest of constraints). In the second order, PSS/Dirac, formulation, 
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canonical transformations also affect the primary constraints; but the form-invariance of the 
whole algebra is preserved (see more details in j^Q, 2] and discussion in next Section). 

In this Section we considered two operations: using different linear combinations of ter- 
tiary constraints and canonical transformations that involve only non-primary variables. 
What is the relationship between these two operations? We showed that they are inde- 
pendent in the following sense. Considering different choices of tertiary constraints for the 
same canonical transformation produces different algebra among constraints, but preserves 
its form-invariance. For example, compare f °° (1841) with r°° for the canonical transformation 
given by ( 1931) . ( 19^1) . ( 1961) . ( |98l) . The PBs for them are different (compare, for example, ( |86l) 
and (189"j) ); but the form-invariance is preserved in accordance with ( I109l) - (lll2p . Applying 
different canonical transformations to the same choice of tertiary constraints modifies the 
constraints and structure functions, but it also preserves the form-invariance of the algebra 
of constraints. We briefly discussed examples of such operations with tertiary first class con- 
straints and canonical transformations that did not involve primary variables (and primary 
first class constraints). We will consider operations with them in the next Section. 



V. THE ROLE OF PRIMARY FIRST CLASS CONSTRAINTS 

Let us discuss the special properties and distinct role of primary first class constraints 
in the Hamiltonian formulation of gauge invariant theories. In the previous Section we 
discussed and demonstrated, by examples, the independence of two operations: the choice of 
combinations of non-primary first class constraints and canonical transformations of phase- 
space variables. We will show that these two operations are not independent any more for 
primary first class constraints, which are either pure canonical momenta as, for example, 
in the Hamiltonians of the affine-metric formulation and the metric formulation due to 
Dirac [3fl , or pure momenta plus some extra contributions, as in the metric formulation due 

n 

to Pirani, Schild and Skinner (PSS) [13( (the oldest one). In the PSS formulation, primary 
constraints are originated from terms in the GR Lagrangian linear in the temporal derivatives 
("velocities") of the go^ components of the metric tensor, and in the Dirac formulation or 
any first order formulation of gauge invariant theories (e.g. affine-metric or tetrad-spin 
connection), from variables without temporal derivatives in the corresponding Lagrangians. 
Momenta conjugate to such variables are primary constraints and at the same time they 
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are part of a phase space of a considered system. This part of a phase space is often and 
mistakenly neglected in Hamiltonian formulations of GR (e.g. see discussion on p. 47 of J] 
and references therein). 

In the Dirac approach to the Hamiltonian formulation of constrained systems, all vari- 
ables are treated on an equal footing and each variable has the corresponding momentum. 
Moreover, variables that are often neglected in the Hamiltonian formulation of GR even 
have a special name given by Bergmann: "primary" , that reflects their importance. In 
monographs on constrained dynamics, e.g. |9|, llOj, and in non-GR Hamiltonians (Maxwell, 



Yang-Mills), primary constraints are always present and are part of the total Hamiltonian 
(name given by Dirac 8|). 

For the Hamiltonian of first order metric-affine GR any change of the primary first class 
constraints would be artificial as they are already in the simplest possible form. So, to 
discuss canonical transformations that involve primary first class constraints we refer to two 



Hamiltonian formulations of the same theory. 



functions of their algebra are quite different 



'SS and Dirac. Constraints and structure 



14j ; but in both cases the complete sets of 



first class constraints lead to the same gauge invariance, as it should be. And this invariance, 
derivable from the constraints, is the four- dimensional diffeomorphism that follows directly 
with no need for a field dependent redefinition of the gauge parameters. For the ADM 
formulation the gauge transformations differ from four-diffeomorphism and the only so-called 
"correspondence" with diffeomorphism 33], or "diffeomorphism- induced" 38| . or "specific 
metric-dependent diffeomorphism" [3^], etc. can be accomplished. 

^he two total Hamiltonians of PSS and Dirac, H^ ss and H® irac respectively, are given 



by 



3 



and 



H? ss = g 0p (K p ss ~ <P° P g)) + g 0p H°pss 9^) > 



113) 



H. 



Dirac 



• Op . n/Op I km „ \ 

90pK Dirac + 90pl^Dirac V 1 > 9pu) ■ 



(114) 



In 15| it was explicitly shown that the phase-space variables of the two formulations are 
related by a canonical transformation that can be performed if one wants to simplify the 
primary first class constraints of the PSS formulation. Note that Dirac [3] found a suitable 
change in the EH Lagrangian by adding to it two total derivatives, which does not affect the 
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equations of motion, but the primary constraints can be brought into simple form. In [15| it 
was shown that the same simplification can be accomplished at the Hamiltonian level using 
canonical transformations in its phase space and, of course, with the same result. What 
is important is that the comparison of two formulations [l5| with different constraints and 
structure functions in the algebras of constraints gives exactly relations (I109l) - (lll2p . Note 
that primary first class constraints of the Dirac and PSS formulations have non-zero PBs 



with a particular choice of secondary constraints that were used in [14| and 4|; but this 
part of the algebra and structure functions also preserves form-invariance under a canonical 
transformation. So, conditions f ll09p -( |TT2|) are satisfied by all first class constraints; and this 



is in complete correspondence with the Dirac conjecture [8| that all first class constraints 
generate gauge symmetry. It is to be expected that all relations amongst them must be 
preserved under canonical transformations to keep invariance in tact. 

As in a first order formulation, different combinations of secondary constraints for the 
Dirac Hamiltonian ( 11141) can also be considered; and one particular combination was dis- 
cussed in {4] 

H T = g o y p + (-g°T 1/2 H- ^fU ■ (H5) 

Note, there are some similarities with our first choice of tertiary constraints for part H" 
( |83|) . One particular choice, out of many possible linear combinations of constraints, as in 
the case of first order formulation (where a few choices of tertiary constraints were con- 
sidered in previous Section), is not special and cannot affect the gauge invariance of this 
Hamiltonian; so it has to lead to a complete restoration of four-dimensional diffeomorphism 
with field independent gauge parameters as was stated in [4j. Of course, the algebra of these 
combinations is different compared with algebra of constraints for (11141) ; but all the different 
choices should not affect the physical results. Similar combinations can also be constructed 
for (11131) ; and under canonical transformations, the form-invariance is also preserved, which 
is the same interplay of linear combinations of non-primary first class constraints and canon- 
ical transformations that we illustrated for tertiary constraints in previous Section. 

Now we briefly discuss a connection of ( I115D with the conventional formulation due to 
Arnowitt, Deser and Misner (ADM) [2j in which the total Hamiltonian is (see e.g. [161 ]) 

H T = NP + N i P i + NH + N i H l (116) 
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where the secondary constraints are exactly the same as Dirac's and are known as the 
"Hamiltonian" 7i and the "spatial diffeomorphism constraint" Tii. According to Pullin 



40| | "It [ADM paper [4l|]] bases the formulation on the Palatini action principle", which is 
actually the affine- metric formulation of Einstein l| (see footnote 3). But it is clear that it 
differs from our results for the affine- metric formulation (I2"9"j) . (1601) and rather obviously have, 
at least, some similarities with Dirac's second order formulation (11151) where the coefficients 
in front of constraints for one out of many possible combinations were redefined and called 
new variables: 

iv=H°y 1/2 , air) 

which are known as "lapse" and "shift" functions 4 . 

As in the examples from the previous Section, because we are working in a phase space 
we cannot just use some invertible transformations like (I117l) - (lll8p and simply write the 
new total Hamiltonian as (11151) . To preserve canonicity, the change of variables (I117I) - (I118I) . 
which is obviously invertible, has to be accompanied by a change of the rest of phase-space 
variables or, at least, some of them (as in the examples considered in previous Section). 
But, as in ADM approach, the space-space components of the metric tensor, gk m , are not 
changed and are exactly the same as in the Dirac formulation 

{9k m ^ km ) Dirac = {g km ^ km ) ADM ; (119) 

and the search for new momenta has to be restricted by 

g Qp K Qp = NP + N*Pi . (120) 



4 Note that even the names used in this formulation manifest their non-covariant nature and shows the 
distinction of these variables, or different roles that they play in the ADM formulation. Their names 
appeared soon after the original works of ADM were published. To the best of our know ledg e, it was 
Wheeler who coined the names of these variables, "lapse" and "shift" in [33J. DeWitt in [29( reserved 
the name "Hamiltonian constraint" only for TL as it is a "particularly important constraint" , probably, to 
reflect its distinction from TLi, although TL and TLi are both the part of the Hamiltonian. 



32 



But it is impossible to find such a transformation that preserves (I120p with the additional 
condition fll 19j) . This unavoidably leads to a conclusion that passing from f 1 1 1 5 j) to fll 16j) 
is not a canonical transformation (an interested reader can find more detail in Section 4 of 
4]). The calculation of one simple PB is enough to prove non-canonicity (see Eq. (152) of 

{N,n km } = j(-3 00 )~ 1/2 ,7r fcm } ^ o. (121) 

So, because the ADM Hamiltonian and Dirac Hamiltonian of GR are not related canoni- 
cally, any connection between them is lost; but the disappearance of four- dimensional diffeo- 
morphism cannot be explained just by this fact. Even if the ADM formulation is considered 
as a model, not related to the Hamiltonian formulation of metric GR, one can argue that 
it still might have diffeomorphism invariance and base his arguments on the fact that the 
number of the primary first class constraints are still the same, and this number defines the 
number of gauge parameters. Note, if the ADM Hamiltonian (11161) is treated as a model, 
then the lapse and shift functions are canonical variables of this formulation. But if one 
claims that (jllfip is a canonical formulation of GR, then it is not the case because ADM 
variables are not related canonically to the metric tensor and its momentum. To accommo- 
date these two possible understandings we will use quotation marks for "canonical" ADM 
Hamiltonian 5 . 

It is not possible to find a canonical transformation, part of which constitutes fll lTf) -( TTT8t) 



with the corresponding momenta (11201) and which simultaneously preserves the condition 
that the space-space components of g^ v remains the same as in the Dirac formulation (11191) . 
But by relaxing this too restrictive condition (I119p . the canonical transformation can be 



found and its form was given in {4] (see Eqs. (156)-(158)) and the problems that arise 
with such transformations were discussed [4|. Note that such a transformation, of course, 



5 The common statements as in 42] "Unfortunately, the canonical treatment breaks the symmetry between 
space and time in general relativity and the resulting algebra of constraints is not the algebra of four 
diffeomorphism" has a double meaning. If "canonical" is understood as the formulation of ADM this is a 
true statement; but for GR, for which ADM formulation is not canonical, Pullin's statement is wrong. We 
would like to note that in the Hamiltonian formulation of GR the covariance is not manifest, however, it is 
not broken as the gauge symmetry of GR, four-diffeomorphism, is recovered in manifestly covariant form 
(see 0, Ql)- So, the Hamiltonian formulation of GR does not break the main property of the Einstein 
GR: general covariance. 
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converts the Dirac Hamiltonian into a form which is different from ADM anyway. So, for 
systems with constraints, even canonical transformations can lead to some problems; and 
for such systems the canonicity of the transformations is the only necessary condition to 
have equivalent formulations. Note that our conclusion based on a particular model, is in 
contradiction with general discussion of canonical transformations for constraints systems of 



431 ] where the authors stated that the condition of canonicity is "too strong" for constraint 
systems; but in our opinion, based on a particular theory, it is too weak, at least for covariant 
theories with first class constraints. 

One obvious, at least for covariant theories, problem with the ADM variables is in the 
original transformation (111 7ft - (1 11 81) and it is related to simple dimensional analysis and to a 
special role of the primary constraints. If we can find canonical transformations, they have 
to preserve f)120p so, in particular, we will have the following part in Ht 

H T = NP + N i P i + ... (122) 

with simple primary first class constraints P and Pi. However, if fields have a physical 
dimension, then the components of the metric tensor should have the same dimension. It is 
obvious that the lapse and shift functions defined by (I117l) - (lll8p have different dimensions in 



terms of the dimension of the metric tensor (dimiV* = and dim N =1/ ydimg 00 ); so the 
corresponding momenta in (11221) (which are primary constraints in such formulation) should 
have different dimensions as well. And such a dimensional mismatch of primary constraints 
for a covariant theory guarantees the failure of this formulation to preserve covariance and 
so such a transformation, even being canonical, has to be rejected in any Hamiltonian 
formulation of a covariant system. This conclusion is based on the following arguments 



related to the special role o 



the Castellani procedure 



primary constraints in the derivation of gauge invariance. In 
16], which we used to restore four-diffeomorphism invariance for 
the PSS and Dirac formulations 14j and 4j, the gauge generator is started from primary 
first class constraints 

G = d^ ) e»P ll + ... (123) 

where is the temporal derivative of n-th order of the gauge parameters n depends 
on the number of generation of the constraints (if secondary constraints are present then 
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n = 1, if tertiary: n = 2, etc.) For a covariant theory should be a true four-vector. 
Again based on dimensional analysis, if the primary constraints have a different dimension, 
then the components of the gauge parameters also have different dimensions. It is the well- 
known fact that four quantities combined together do not necessary form a true four-vector 
as their components must transform in the same way under general coordinate transforma- 
tions 44j. As we can see from (11171) and (I118p . iV and N l transform differently, because 
they are defined in terms of particular components of the metric tensor in non-covariant 
way (see again (I117P and (I118P ). This property is also transferred to the corresponding 
momenta and gauge parameters. If, at least one gauge parameter has a different dimension 
or transforming property from the remaining parameters (which is exactly the case here), 
it is just impossible to combine them into the four-vector gauge parameter, which is needed 
for four-diffeomorphism. To conclude: even in the case of canonical change of variables, but 
with a mismatch of their dimensions (at least for variables which correspond to primary 
constraints) the covariance is lost. So, introduction of lapse and shifts functions by itself, 
whether they are a part of a canonical transformation or not, unavoidably destroys covari- 
ance and in turn the equivalence with the original covariant theory. Please note that there 
are different approaches to the restoration of gauge invariance, where a generator is built 
on other principles and is actually started from the end of the constraint chains (e.g. see 
[12)]). i.e. from secondary constraints (in second order formulation) multiplied by the gauge 
parameters 

G = ... + eft + eUi . (124) 

In this case the conclusion is the same, because if lapse and shift functions have different 
dimensions, so do the "Hamiltonian" and "spatial diffeomorphism" constraints and the cor- 
responding gauge parameters in (j!24p . We will discuss in detail the methods of restoration 
of gauge invariance in Sections 6 and 7 with application to a first order affine-metric formu- 
lation of GR. Here we would like just to add, that in both methods, the initial assumption 
of either the Castellani algorithm or HTZ ansatz is the independence of gauge parameters of 
fields; and this assumption is used in the iterative procedure to a find generator. So any field 



dependent redefinition of gauge parameters, as advocated in many articles (i.e. [3 71. l38l. l39|). 
that is performed after completion of the procedure, is in complete contradiction with the 
initial assumptions for both approaches. 
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Now we will consider one additional example, the Hamiltonian formulation of the 
Einstein-Cartan (EC) theory, which is a little bit aside of the main topic of this article. 
But it provides an illustration of even further (possibly general) restriction on the manip- 
ulations with primary constraints (canonical variables), which cannot be illustrated using 
metric or affine-metric formulations. 

The Hamiltonian formulation of the Einstein-Cartan theory in its first order form, the 
so-called tetrad-spin connection formulation was discussed in many works (i.e. Q|, 4?|] ) . 
As in the Hamiltonian formulation of affine-metric GR, it leads to second class constraints 



that should be eliminated, and after the Hami 
Hamiltonian (up to a total spatial derivative) 



tonian reduction leads to the following total 



45 



4g| 



H T = e 0{p) n ^ + u (ap)7i oia(3) + e mX 0{p) + Uo {a 0) X° (a(3) ■ (125) 
In the 3D case, to obtain (j!25p is a simple task because there are no secondary second class 



constraints 



48| . In the AD case, the different methods (specific to this dimension) of solving 
secondary constraints were used. One particularly transparent method is the introduction of 
Darboux coordinates (specifically constructed only for AD) due to Banados and Contreras 



49( | . Of course, the Dirac procedure can be used in any dimension higher than two, 6 which 



gives the same Hamiltonian (11251) . and that was shown in 45|. Direct calculations are 



involved and a considerable simplification occurs when the Darboux coordinates (common 

n 

to all dimensions) are used [501 ] . 

The notation can vary from paper to paper, but it is usually explained in detail. The 
equation (11251) is written in the notation used in Q|; and it is very close to what can 
be found in the first "gauge-free" formulation 7 of the Einstein-Cartan Hamiltonian, due to 

n 

Castellani, van Nieuwenhuizen and Pilati [46], where the canonical variables (after elimina- 
tion of secondary second class constraints) are 

e^ p) ,7i^\u 0{af3) ,7i ^l (126) 



6 As for metric EH action, the second and the first order, affine-metric, formulations are equivalent only 
in dimensions higher than two, the same thing happens for EC action: tetrad and tetrad-spin connection 



formulations are equivalent also when D > 2 (see 45j . Section II). 



7 "Gauge-free" means without fixing a gauge at the beginning of analysis. Such a fixing is in contradiction 
with the Dirac procedure as a gauge cannot be fixed before a gauge symmetry is found. 
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The form of (I125p . which is a nice covariant expression for the total Hamiltonian of the 



first order EC action, was known 
with a few rare exceptions (e.g. 



or a long time. But it is difficult to find it in this form, 



5l| where the first two terms of (I125p are given in Eq. 



(2.4) and the last two in Eq. (3.3)). Because the "canonical" formulation, in accordance 
with the conventional wisdom, is equivalent with the presence of lapse and shift functions, 
the change of variables is always performed to introduce them 8 instead of completion of the 
Dirac analysis (proof of closure) and restoration of gauge invariance for (11251) . (Some steps 



in the analysis of (11251) can be found in 



45 



■) 



The Hamiltonian in terms of lapse and shift functions for tetrads becomes (this form is 
much easier to find in literature, contrary to (I125p . e.g. [521 ]) 



H T = NP + N l P t + u 0[aP) ^ + NH + mU + u 0{aP) x 



Ma/3) 



(127) 



Note that introduction of lapse and shift functions guarantees the disappearance of co- 
variance in the formulation. We have already discussed the effect of the ADM variables, 
which have different dimensions; and any hope to have a covariant formulation is lost. They 
also, as in metric formulations, are not canonical for a transformation from (I125p to (I127p 
if they are introduced with restriction on the rest of tetrad components and corresponding 
momenta (which is always the case). 



The canonicity of transformation in the phase space, i.e. 



rom the complete set of 



canonical variables of (TT251) (e^ p ), ^ {p \ uj (a/3), vr 0(a/3) ) |45|, |46|, [511] to variables of (1127) 
(N, N\ P, Pi, efc(p), n k ( p \ ^o(a/3)> vr 0(a/3 )), has never been discussed and it has the same defi- 
ciency as for passing from the Dirac to ADM formulations in the second order metric EH 
action. Moreover, exactly as in the metric case, only the part of variables is involved in such 
a change 



eo(p),7r 



N,N\P,P { , 



(128) 



(e fc ( 



IT 



Hp) 



^0(0/3); 



7T 



0(a/3)^ 



)EC = {e k[p) y {p \uj^ ) y^) 



ADM 



(129) 



where the lapse (11171) and shift f II 1 8 j) functions have to be expressed in terms of the original 



Such a transition is justified by either "it is more convenient" 51] or "it is useful" 4 
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phase-space variables (I126p of the reduced Hamiltonian (j!25p using g^ u = e'f ■ ) e u ^ 



^=H?P) e0(p) )" 1/2 ' ( 13 °) 



N l = — 7r~~nH\- (131) 

e O pf o iP ) 



'(-rV 

If 

As in the metric case, there are no canonical transformations for the subset of phase-space 



variables (I128p if the rest of variables of (11261) is not involved. As in the transition from the 
Dirac to ADM variables it was enough to find one PB (11211) which proves non-canonicity of 
the ADM variables, for the tetrad formulation at least one PB is also non-zero 



(132) 

So, introduction of lapse and shift functions in the tetrad formulation is also a non-canonical 
transformation. Actually, based on previous analysis, it is obvious that such a formulation 
(even with adjustments for canonicity that will change constraints) will create a dimensional 
mismatch of variables and gauge parameters, and so destroy the covariance in exactly the 
same way as lapse and shift functions destroy it for the metric GR. The combination "Dirac- 
ADM" is not correct, as well as "ADM-Einstein-Hilbert" , so the ADM and Einstein-Cartan 
formulations are not compatible. All problems related to the "spatial diffeomorphism" 
constraint safely propagate into the tetrad formulation of GR and such the Hamiltonian ( 11271) 
is not the Hamiltonian of the original theory. In addition, this new theory is not covariant 
by construction. However, the use of ADM variables for first order tetrad-spin connection 
formulation is much more interesting example compared to its metric counterpart; and this 
is the main reason to include the Einstein-Cartan Hamiltonian in our discussion about the 
role of primary first class constraints. 

Here one can make an additional and simple observation related to the role of primary 
first class constraints. Constructing the generator for formulation ( I125p . one obtains 

G = e {p) n°^+ 6^)71°^ + ... (133) 

with two gauge parameters, "rotational" (£( a m) and "translational" (s( p )), both with "inter- 
nal" indices (which correspond to motion in the tangent space], and which lead to transla- 
tional and rotational invariance in the internal space in 3D 48J, [53(, as well as in all higher 
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dimensions [45l . [5J| . It must be emphasized that among the many papers on the Hamilto- 



nian formulation of EC action, the work of 



461 ] is an exception because this is the only one 



where lapse and shift functions are just a short-hand notation, so non-canonical transforma- 
tions f)130p . (11311) were not performed. Should one apply the Castellani procedure (see next 
Section) starting from the primary constraints of |46j the gauge invariance, translation and 
rotation in the internal space, would follow (not a diffeomorphism, which is a symmetry, but 
not the gauge symmetry of the EC action). 

In contrast, if the generator is built for (11271) 52J (let us forget for a moment about 
non-canonicity) then it becomes 



G = eP + e'Pi + e {a ^7T ^ + ... (134) 

which also has two gauge parameters. One, as before, corresponds to rotation in the internal 
space (£( a /3)); but the second (actually, there are two of them, e and e l ) corresponds to the 
translation in the external space: e for time lapse and e l for shift in a space-like surface. 
Because P and Pi are momenta conjugate to lapse and shift functions, which have different 
dimensions and rules of transformation, then P and Pj also do not form a true four-vector, 
and so are the parameters e and e l (it follows from (11341) ). Even if by relaxing the conditions 
of (11291) and finding some canonical transformations, and even with a match of dimensions, 
such a change of gauge symmetry, from translation in an internal space to lapse and shift in 
an external space would be strange. 

Is it possible to have equivalence in such a case? This would mean that the Hamiltonian 
formulation does not give a unique gauge invariance; but this conclusion is hard to accept. 
Why does one, having the simplest possible constraints presented in dl 25j) . have to perform 
some manipulations, besides the desire to have a "canonical" formulation, i.e. to have lapse 
and shift functions that destroy covariance? We think that for field theories, changes of 
variables (even canonical, if we can find such) should be restricted by the requirement to 
also preserve a tensorial character of primary variables. Maybe, such a condition is equivalent 
to the preservation of form-invariance of the algebra of constraints (I109l) - (lll2p . There are 
still many questions that need to be clarified. Dirac introduced his procedure as an outline 
of a general approach, together with his famous conjecture about the connection of first class 
constraints and gauge invariance. Only after the methods of restoration of gauge symmetry 
were developed, it became clear that the Dirac conjecture is correct and it can be used as a 
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procedure. But there are still many other questions, especially related to field theories, that 
remain to be answered. 

Based on the examples considered in this Section, we can make a conclusion that for 
constrained field theories the canonicity of change of phase-space variables is a necessary, 
but not a sufficient condition. To keep an equivalence of two Hamiltonian formulations, the 
canonical transformations that lead to the mismatch of the dimensions for primary variables 
(which is permissible for mechanical systems or field theories without first class constraints); 
and also for canonical transformations that change tensorial character of primary variables, 
must be disregarded. Primary constraints are the true Masters of the Hamiltonian formu- 
lations of gauge theories. They do not need a master constraint programme as they are 
primary constraints; and manipulations with them are restrictive, especially for covariant 
theories. In the formulations considered here where primary constraints are in the simplest 
possible form (e.g. (rQ ,n° )), it is better not to change them at all, to keep in tact the 
primary properties of the system (or there should be a very good reason to do the changes). 

The reasons behind converting covariant expressions for the total Hamiltonian of the 
second order metric GR H 11 3 f) - ( TTX5T) and the first order tetrad-spin connection formulation 
(11251) into non-covariant expressions (11161) . and (11271) . which destroy any hope of having 
covariant results (e.g. four-diffeomorphism for metric GR), remains a mystery to us. (We 
prefer to avoid speculation and leave it for the History of Science to figure out what happened 
50 years ago, and why Einstein's general covariance became unimportant.) But what is even 
more mysterious is the unquestionable acceptance of the ADM formulation by the majority of 
practitioners despite its inconsistencies and despite (not often) the strong voices that express 
concern about its deficiency, contradiction with GR, and even give "hints" about the source 
of the problems. We provide just a few such statements. Hawking 30 years ago concluded 



55( | : "The split into three spatial dimensions and one time dimension seems to be cont rary 



to the whole spirit of relativity." There appear more recently the statements of Pons 39] : 
"Being non-intrinsic, the 3+1 decomposition is somewhat at odds with a generally covariant 
formalism, and difficulties arise for this reason" and Rovelli 56(: "The very foundation 



of general covariant physics is the idea that the notion of a simultaneity surface over the 



ems m 



universe is devoid of physical meaning". The warning about the sources of prob 
further constructions based on the ADM formulation was given by Landsman j|37]: "the 
lack of covariance of the ADM approach, which is especially dangerous in connection with 
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quantum field theory" . 



58j "Thus the full group of 



Finally, we repeat the "hint" given by Isham and Kuchar 
spacetime diffeomorphism has somehow got lost in making the transition from the Hilbert 
action to the Dirac-ADM action" . This statement (without long calculation or a long chain 
of logical constructions) immediately leads to the conclusion that if a transition from one 
action (Einstein-Hilbert) to another (ADM) was performed by a change of variables and 
something got lost then "somehow" is exactly the change of variables that were performed 
in such a transition. We strongly object to the use of combination "Dirac-ADM", as nothing 
"got lost" in the Dirac Hamiltonian formulation based on the Einstein-Hilbert action j^]. 

Last year, 50 years after publication of the Dirac Hamiltonian of GR {3J] 9 , it was ex- 
plicitly demonstrated 4j that his Hamiltonian of the second order metric GR leads directly 
to four-diffeomorphism without any, even numerical, redefinition of gauge parameters. This 
clearly demonstrates that four-diffeomorphism has not "got lost" if one abandons the idea of 
making a transition from the Einstein-Hilbert action to ADM action (i.e. if the ADM vari- 



ables "got lost" instead). In J] we argued (see Eq. (163) of [4j) that, based on equivalence of 
Hamiltonian and Lagrangian methods, it should be possible to demonstrate non-equivalence 
of the ADM and EH actions also at the pure Lagrangian level (without going to a phase 
space); but we did not formulate the criteria and did not show a proof. Recently we demon- 



strated that the Einstein-Cartan action is invariant under translation in a tangent space 



54) 



using the pure Lagrangian method (to verify the strong indication of the presence of this 



gauge symmetry in the Hamiltonian formulation [451] ). This also allows us to formulate the 
condition for equivalence of two actions at the pure Lagrangian level for singular systems: 
if transition (field redefinition) from one singular action (e.g. Einstein-Hilbert) to another 
singular action (e.g. ADM-inspired tetrad) is performed by an invertible change of variables 
(the necessary condition for equivalence) we can find differential identities for both of them 
and obtain the corresponding transformations for their sets of variables. These transforma- 
tions, for equivalent formulations, must be derivable one from another by using the same 
(or inverse) redefinition of fields as in transition from one action to another 10 . 



9 We are thankful to Shestakova [59( for this observation and in this article we are also trying do keep track 

of important anniversaries. 
10 The "loss" of four-diffeomorphism in the ADM formulation was demonstrated by Banerjee at al [60( 
(see Section 6 of [s^]) where symmetries of the ADM Lagrangian and the Einstein-Hilbert Lagrangian 
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This year, one year after 50 anniversary of Dirac's paper j^], it was a celebration dedicated 
to another paper on the same subject |2j that took place in Texas A&M University on 
November 7 and 8, 2009 "ADM-50: A celebration of current GR Innovation". This event 
marked 50 years of loss of four-diffeomorphism and covariance in the "canonical" formulation 
of covariant GR. In the above mentioned article [2| the authors claimed to consider the 
Hamiltonian formulation of GR in Palatini form, i.e. the formulation of Einstein This 
is also the main subject of our paper, but with different results and without loss of four- 
diffeomorphism (see next Section). 

The transition from the EH action to the ADM action, which is responsible for a loss 
of four-diffeomorphism, continues to propagate into new fields. A particular example is 
Loop Quantum Gravity (LQG) - one of the major players in the quest for quantization of 



gravity. According to Thiemann 



621 ] "One of the reasons why LQG is gaining in its degree 



of popularity as compared to string theory is that LQG has 'put its cards on the table'". 



Let us look at the LQG "cards" that were put together in the recent article [63j (see first 
paragraph of Introduction) that describes the meaning of LQG: "...LQG is a mathematically 
rigorous 11 quantization of general relativity (GR)... It is inspired by the formulation of GR 
as a canonical dynamical theory... The total Hamiltonian of GR is a linear combination 
of the Gauss constraint, the spatial diffeomorphism constraints and the Hamiltonian con- 
straint. Thus the dynamics of GR are essentially the gauge transformations generated by 
constraints" . Is it canonical formulation of the covariant GR if the "spatial diffeomorphism" 
constraint is present? It is (11271) that cannot have covariance and it is not equivalent to 
the Einstein- Cart an GR. The Hamiltonian of LQG can be written in different variables; 
but they are originated from lapse and shift functions, that is why, the "spatial diffeomor- 



were considered. It is not a surprise (taking into account equivalence of Lagrangian and Hamiltonian 
methods) that the same gauge invariance as in the case of the ADM Hamiltonian was found also for 
the ADM Lagrangian. The transformations found using the ADM Lagrangian are not convertible (using 
the original change of variables) into four-diffeomorphism invariance of the EH Lagrangian obtained by 
the same method 61|. Only after the same manipulations with a field dependent redefinition of gauge 
parameters 37] (which should also be field independent in construction of the generator from differential 
identities) invariance of the ADM Lagrangian can lead to "correspondence" with diffeomorphism. The 
authors of 6fJ call this "the equivalence between the gauge and diff parameters by devising of the one to 



one mapping . 

11 In our taste, "rigorous" in combination with mathematics is a tautology, but because the rival of LQG 
employs elegant mathematics, the use of just mathematics, probably, sounds a little bit weak. 
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phism" constraint is always there, as in the ADM gravity (11161) . Moreover, and this was 



already demonstrated (e.g. 64j, |65|, [66]), all new variables which are used in LQG can be 



converted into ADM variables by canonical transformations, i.e. new formulations are equiv- 
alent to the ADM formulation. The gauge transformation generated by constraints of the 
ADM-inspired Hamiltonian of tetrad gravity is the "spatial diffeomor phi sm" ; but the trans- 
formations generated by the Hamiltonian of EC action are different I45I . [54J] . So, in reality 
LQG, using rigorous mathematics, studies quantization of the ADM-inspired model and its 
gauge invariance and, because it is not covariant, all well-known results follow: quantiza- 
tion of three-dimensional surfaces |67|. Lorentz violation, etc. The suspicion of people not 
working in this field about how a covariant theory of GR can lead to such results is correct. 
The explanation of these effects originated from the classical Hamiltonian mechanics and 
the theory of canonical transformations. These are no quantum effects, but they just man- 
ifest the dependence of quantum effect on a classical background, the ADM Hamiltonian. 
To show this, neither rigorous nor elegant, but just a mathematical condition, the Poisson 
brackets (11211) and (11321) . that must be zero to keep two Hamiltonians equivalent and the 
notion of a phase space are needed. This is what any outsider to this field can and should 
read from the LQG "cards". This situation can be perfectly described using one of the 
famous paradoxical forms of Wheeler (e.g. see 68(): LQG is quantum "gravitation without 
gravitation" in Einstein's sense. 

We discuss recent, current and coming soon (see next paragraph) anniver- 
saries/celebrations. But let us also mention a lesser known anniversary: 200 years ago 
Simeon-Denis Poisson published his work "Sur la Variation des Constantes arvitraires dans 
les questions de Mecanique" in the Journal de l'Ecole Polytechnique, Tome VIII, p. 266, 
Decembre 1809, Paris. This is the first appearance of what has become known as the Poisson 
brackets; and we have used them to show that transition from EH or EC variables to ADM 
variables is not canonical (see ( 11211) and ( 11321) ). 

Let us return from the past to present days, and continue with one more quotation from 
Thiemann 62J: "The 'rules of the game' have been written and are not tinkered with". And 
that is an important standard to which current research in GR must be held. The rules of 
the game are those published by Poisson, in 1809, and by Einstein at the beginning of the 
twentieth century. These are important rules which must not be tinkered with, and must 
not be ignored. 
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However, seems to us that the main rule of the LQG "game" is to work with the ADM- 
inspired action which is not covariant and, because of this, allows to obtain non-covariant 
results but, at the same time, to present this as properties of the Einstein-Cartan action, 
as scientific community at large still associates gravity with the name of Einstein. Neither 



outside view on LQG (Nicolai, at al [69]) nor inside view on the same subject (Thiemann 
6^ |) disobey this rule, and, as a result, LQG is also approaching its 25 anniversary. 

Of course, there is a freedom of research (or variety of "games" in nowadays language) 
and one can choose to play with the ADM-inspired metric or tetrad models in a framework 
of non-covariant theories; but one should not be surprised by the results which are not 
covariant. Or one can study the covariant theories of GR (Einstein or Einstein-Cartan) and 
in this case one ought be surprised if non-covariant results follow. 

We choose to work with the covariant Einstein formulation of GR and do not convert his 
brilliant theory into "common currency" 12 , i.e. we do not make any non-canonical change 
of variables and we return to the total Hamiltonian ( |29|) . ( |60l) of the first order affine- 
metric Einstein GR, the main subject of this article. The Hamiltonian formulation of any 
system leading to first class constraints cannot be considered as complete without deriving 
the corresponding gauge invariance. In addition, a restoration of gauge invariance is an 
important consistency check, especially if we investigate such a complicated theory as GR. 

In the course of our calculations it becomes clear that, despite different starting points 
in our current consideration (affine-metric action (T5J)) and those based on different set of 
variables (see (ED) that were discussed in 0,0 , and especially due to work of Ghalati 
and McKeon 7] where the closure of Dirac procedure was demonstrated for the first time, 
the results are very similar. One can say it is remarkable or surprising, but nothing here is 
either "remarkable" or "surprising" because these results come from two equivalent first order 
formulations of the same theory (see Appendix of [5j). Of course, there are some differences; 
but secondary and tertiary constraints (for the same particular choice) are exactly the same 
after the canonical transformations (|93l) -( l98l) were performed. Neither our formulation, nor 



[ll were converted into a "canonical" form (despite that some attempts to make a formulation 
"canonical" were made in the novel approach of [ill); and one should expect the complete 



12 The name given by Pullin 40] to the ADM formulation in his Editorial note for republication of [4l| in 
70|. 
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restoration o 
metric GR 



four-diffeomorphism as it was obtained for the second order formulation of 



The first steps of the restoration of gauge symmetry were made in llj, with a truly 
surprising result for a covariant theory: the need for a field dependent redefinition of gauge 
parameters to make the transformations of h 00 (the only one that was calculated 13 ) "corre- 



spond to diffeomorphism invariance" 
lations of metric and tetrad gravities 



, exactly as in the conventional "canonical" formu- 



37 



38. 



39]. Why does the Hamiltonian formulation 



of the first order EH action, obtained without non-canonical changes of variables, give so 
different result compared to the Hamiltonian formulation of the second order of GR? In the 
restoration of gauge invariance, we used the Castellani algorithm 16j but in llj the different 
method, the HTZ ansatz [jjj], was employed. This is the only difference, and a comparison 
of the two different methods is the best point to start the search for understanding of the 
apparently different gauge transformations for first and second order formulations of GR. 
The field dependent redefinition of gauge parameters of these two formulations could be just 
an artifact of a particular method (if in one it is needed, but in another it is not). In the 
next Section, we first apply the Castellani algorithm to the first order formulation to obtain, 
as in [ll|] , partial transformations (which is enough to make a conclusion about the necessity 
of a field dependent redefinition of gauge parameters). And then, we compare this with the 
results obtained using the HTZ ansatz. 



VI. CASTELLANI ALGORITHM 

In his paper [lfl Castellani illustrated the application of his method by considering Yang- 
Mills theory (the system with only primary and secondary first class constraints). We 
also used this method to restore gauge invariance in the Hamiltonian formulation of the 
second order metric GR QQ and of the first o r de r , t e tr ad- Spi „ connection, foliation of 
Einstein-Cartan action in the three dimensional case k8j. In the Hamiltonian of the first- 

n 

order affine-metric GR, we have tertiary constraints; but the procedure of [16[ is general 
and can be applied to systems with any number of generations of constraints. We are not 



13 To avoid any confusion, we would like to emphasize that converting the results of direct calculations 
into the ADM form, which is extensively discussed in [llj], was not used in the calculation of these 
transformations. 
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aware of an application of the Castellani procedure to a realistic Hamiltonian with tertiary 
constraints, and such an application is interesting by itself. Our main goal in this, and in 
the following Section, is to analyze the appearance of a field dependent redefinition of gauge 



parameters, not a complete restoration of gauge invariance. So as in 



111 ], we will calculate 



transformations only partially but for all fields, to see whether this method produces correct 
terms in the transformations of all fields, contrary to where the transformation of only 
one field was found. 

In the Castellani algorithm [uj the generator of the gauge transformations for the Hamil- 
tonian with first class constraints for the system with tertiary constraints is given by 



G = e"G (0)A1 + e"G (1) „ + e"G (2)A1 (135) 

where are the gauge parameters and e M , are their temporal derivatives. The number 
of gauge parameters and their tensorial dimension are uniquely defined by primary first 
class constraints, so for the formulation considered, the number of parameters is equal to 
the dimension of spacetime, D. The functions Gu^ are defined by the following iterative 



procedure (see Eq. (16b) and for more details see also Section 5 of [16] ]) 

G ( a> = nj°, (136) 

G (1)M + {G (2)M , H T ) = J d^a; (1?, V) < (V) , (137) 

G (0> + {G (1)M , H T } = j d^p; V) n°° (V) , (138) 

{G (0)m , H T } = primary. (139) 

Note that only primary constraints enter equations (I136I) - (I139I) explicitly. The functions 
G( 2 ) M are uniquely defined as primary constraints. The functions, G(i) M and G( ) M , in general, 
are not just secondary or tertiary constraints because, for example, different linear combi- 
nations of tertiary constraints can be considered (see (1531 . (1551 . and (1571) ). This makes this 
method insensitive to our choice of combinations of non-primary constraints and gives the 
same gauge invariance regardless of what combinations we will call 'tertiary constraints'. We 
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will illustrate this important point in more detail: for complicated theories the possibility 
of working with different combinations of constraints, without destroying its unique gauge 
symmetry could give significant computational advantages. 

Using PBs among the first class constraints and the total Hamiltonian, which is given by 



and Gmp becomes a secondary plus linear combination of primary constraints with 
coefficient-functions a" that have to be found. To shorten the notation, we will not write 
integrals in equations that involve coefficient functions, which in general might also depend 
on fields and their derivatives (see, i.e. 4], Q]). Only in the case of contributions with deriva- 
tives, which is specific to field theories, a more careful treatment is needed for finding the 
expressions of the corresponding coefficient-functions. In this Section we restrict ourselves 
to some simple steps of procedure mainly to discuss its important properties and to show 
that it provides a strong indication for the correct restoration of the gauge transformations 
for all fields. We will consider only a part of the transformations that are produced by 
simple contributions from coefficient-functions, which do not involve spatial derivatives, so 
our short notation without integrals will not lead to any ambiguity. The complete restora- 
tion of diffeomorphism invariance from the constraint structure of the Hamiltonian of first 
order affine-metric GR using the Castellani method is in progress and, of course, the com- 
plete calculations involve the spatial derivatives of the fields and integral form of expressions 
with coefficient-functions becomes important. Details of such calculations will be reported 
elsewhere. 

At the next step of the procedure, using (11381) and the above result (11401) the functions 
(j(o)/j can be found. The total Hamiltonian (in our condensed notation) is 



In such a form, fl50l) - fl52l) . H' c is independent of "primary fields" Tq (fields for which 
conjugate momenta are primary constraints) and this allows us to perform a few steps of 
calculation, which are independent of a particular choice of tertiary constraints. Using this 



(Egjl - Pty with (HHD and (j5D} - (j52jl . we can solve (jlHBjl for G (1)fl 




Ht = I4HJ° - TUl + K ■ 




47 



form of the Hamiltonian, (11411) and after some rearrangement, (11381) becomes 



Gm, = - M°, Koxl } + { x ™, h' c ] - n? K, n»} i^ - + n? K, Ho 

-n°°K,^} + ^ n ?°- ( 142 ) 

At this stage of the calculation, both coefficient-functions a 1 ^ and {F enter (1142ft while at 
the previous step only was present in (I140p . 

The last equation of Castellani algorithm, fl 1 3 9 j) . serves to find unspecified coefficient- 
functions oi v and fi v as all terms proportional to the secondary and tertiary constraints 
must be identically zero 

{C? (0)A4 , r^ o n°° - r^ oX °° + H' c } = primary. (143) 
After a v and (3" are found, we have the following generator 

G = e"G (0 )„ + e"G (1)/1 + e"G (2)/i = (144) 

£ ^- jX^ ,X„ I i oo + ,ti c ) - \ a ^ l u ^oo li -a ) + %Xv +P;ll 7 j+e^ y-Xn + j+ e ^ 

that allows us to find the transformations of fields or combinations of fields, F, using 

5F={G,F}. (145) 

Here, we again do not specify an H' c , which can be written in different ways, as its form de- 
pends on a choice of tertiary constraints. Obviously a particular choice of tertiary constraints 
can only affect the transformations through the coefficient-functions. The transformations 
generated by the part of the generator (11441) . without coefficient-functions is 

*F = (- M°> X°°} r ' + {xl°, H>}) - e»xT + ^n? , F} (146) 

and it produces the same result, whatever combination of tertiary constraints is considered. 
The complete restoration of the gauge generator, and especially the gauge transformations 
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of all fields, is a technically involved problem. Solutions of second class constraints must be 
used or we have to go to the reduced Lagrangian fl90l) . which corresponds to the reduced 
Hamiltonian, and find the momenta in terms of the coordinates from the equations of motion, 
as it is described in the HTZ paper jl2j |. We restrict our discussion to the first and relatively 
simple steps of calculations: to single out a method of restoration that does not lead to 
contradictions. 

Let us start from the derivation of the partial transformations of Tqq to illustrate that 
the Castellani algorithm is independent of a choice of tertiary constraints. We will show 
that it leads to the correct contributions to the gauge transformations for all fields; and that 
these transformations are equivalent to diffeomorphism without any field dependent or even 
numerical redefinitions of gauge parameters. 

First of all, we write G(o)/j in components that allow us to explicitly calculate some PBs 

n - r fc v»» _l /v"0 u'X rroo r7 

Lr(0)0 - i QoXk + \Xo > n cf ~ ll u TFT 1 00 

01 00 

- «oX°° + n°° K, x? } - n°° k, h' c } + (147) 

<^{o) P — ~X P 1 oo + \X P ) n c) ~ LL u TFT 1 oo 

01 oo 

- <x°° + n°° { a ;, x ° } rj, - n°° {<*;, h[) + p;n™. (us) 

The condition (11431) gives for fll47l) - fll48p the following terms, which are not proportional 
to primary constraints, 

{G m , H T ) = xTt k 00 - r k ooX ™r 00 + T k 00 {xl°, H' c ) (149) 



- { {x°o°, K) , x° } ^o + (W, K) , H' c ) - ^lt« oX ™ 

01 00 



00 SOq ■ , 00 QOlpa n v f _00 tti\ ^00 f„u tt/\ 

-Xu TFcT oo + \ a oX,y , Xa ) 1 oo ~ a o \Xy , U c } - X v {ot Q , UJ 



^oo 



+ K, x° a } r? x°° - {<, H'} X °° + & 7 x 7 = o, 
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(150) 



-x?°r° + {- x ™r° 0) -r^ } - r°, K) - , xl } n + {{x™ K) , h' c ) 



. 5a v 

_i OOTp^Xi, - Xz/ t a p' n a / 1 00 + \ a pXu iXu J 1 00 ~ a p \Xu > -"c / ~ X* ("p, ^c} 
01 00 



+ K x? } r?ox°° - K, x°° + /^x°° = o. 

The only terms in ( I149l) - fll50p that can give contributions proportional to tertiary con- 
straints (hidden in H' c ) are: 

+ Too W, K) - ^} , x° } r« + { {x™, ^} , ^} - ag {x°°, ^} = 0, (151) 

- rgo {xl?, K) - { «°, , xl } K + {{x°p°, K} , h' c ) - a ; {xl°, K} = o. (152) 

Equations (I151l) - (ll52p allow us to find the coefficient-functions a". Note that ()15ip - (jl52p 
are algebraic equations which are written in a general form and are independent of a partic- 
ular choice of tertiary constraints. But to find the coefficient-functions we have to expand 
the expressions (115ip - (1152p and collect terms proportional to the secondary and tertiary 
constraints. To do this we have to specify our choice. A few choices were discussed in Sec- 
tion 4; but let us start from one particular combination, fl60|) with the simplest PBs among 
constraints, namely one that leads to zero PBs among secondary with tertiary constraints: 
{xl°, T °°} = 0. The simplicity of the PBs was also the reason in [7] to use these combina- 



tions to prove c 



osure of the Dirac procedure. It was precisely this choice that led to the 



conclusion in ll| that diffeomorphism invariance does not follow directly and a field depen- 
dent redefinition of parameters is needed to find the "correspondence" with diffeomorphism 
using HTZ ansatz [12|. So, we specify H' c (note that until this moment all equations were 
independent of a choice of tertiary constraints) and use the following form 

50 



H 'c = 7^0° + + A°x o° + A m xZ (153) 

With this H' c , and keeping only terms that lead to tertiary constraints (needed to find a") 
from (I15ip and f l 1 531) . we have 

9F fc oo , ( J_ r oo , ^. r oo^\ r o _ o ( J_ T oo , h m 7 .oo\ _ n k oo 
Zi oo r fc + \h 0Q h 00 m J 00 V h 00 /i 00 m / * 

+ {^o°° + + M°, A *" + ^ffl , = 0. (154) 

Let us, for simplicity, restrict our calculations even further and consider only the depen- 
dence of the coefficient-functions a 1 ^ on Tq . From the first line of f l 1 541) we can uniquely find 
how «q and depend on Tq (there are no contributions proportional to Tq in the second 
line of (11541) ). Combining together the terms in the first line of (11541) that are proportional 
to the tertiary constraints we obtain 

2r 00 - a + — r 00 - — a r fc + — r oo - — a r =0. (155) 



froo uu h 00 7 V ft. 00 /i 00 

The second bracket of fl!55[> gives 



a °(r) = r° . (156) 
Using this result and the first bracket of (11551) . we obtain 

«o(r) = 2r fc . (157) 

Similarly, from (11521) we find 

o£(r)=0, a k p (T) = -5 k p T 00 . (158) 

With these results we can also calculate the contributions to the coefficient-function fi v 
from the corresponding expressions that are proportional to the secondary constraints. In 
fll49p - fll50p . keeping only terms with fl v and those with temporal derivatives of Tq , we have 

OOf k ™ $ a 00 00 ^ a fa , fl0_.00 , tfj> 00 _ n /-i cqN 

Xfc 1 oo _ 1 ooTF^X,, - Xv TF^r v oo + PoXo + PoXp - U, (lt>yj 

01 00 01 00 
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- x°A - ^K xl° - x°°T^r « + pX° + W = o, (ieo) 



5a u ■ 5a u ■ 

Z_2E_r a -v 00 _ -v 00 p t 

lp A 00 rp« 1 OOAv rp a 1 00 ^ Vp, 

01 00 01 00 

which after substitution of (I156l) - (ll58p gives 

/3 ° (f ) = 2T° 00 , f3 k (r) = 3r o fe , (161) 

/3 p ° (r) = 0, (3 k p (r) = -5 k p t° 00 . (162) 
Let us check the consistency of these simple partial results and consider variation of 5Tq 

^ r oo = { G > r oo) • (163) 

For the relevant part of the generator (terms proportional to primary constraints n°° and 
found in (I156I) - (I158I) and (I16ip - (I162I) Independent parts of and /3^), we have 



g = e° (r°X° + f M° j + e°K n™ + £° 2 rM - £ p KoK° + + ^ P K°- ( 164 ) 

This part of the generator gives 

Ko = ~e°t ° - e°r° 00 - e°, (165) 

ST p 00 = -e°tl - 2i°Tl + e*T° 00 - & . (166) 

We have to compare these partial results with the well-known transformations of the 

components, Tqq, under diffeomorphism invariance. Using the general expression for the 
transformations of affine-connections under diffeomorphism, 

c "pA A I "pp A c-PT^^ c-P c-P ( 1 

°diff L ~ ~ 6 ,pv ' 1 pis 6 ,p ~ 6 1 iw,p ~~ 1 n(f- ,v ~ 1 up 6 ,p ) {XOI ) 

for two particular components of (11651) and (11661) . we obtain 

ddiffT qo = —£°oo ~ Too^o^ — £ °do^oo + Fo d k e° — e k d k To — 2 (S fc + ^Tm) 9oe k , (168) 
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— e°d T b 00 — E k dkT b 00 — 2T b 00 doe° — 2T, 0k d e + 



2 




(169) 



D - 1 



Here we used the redefinitions (@|, (TT5l) and (TT8T) (from last two it follows that ir^ = E M!/ ). 
We do not perform such substitutions for and as they do not affect the parts of the 
transformations that we shall compare with our partial result (I165l) - (ll66p (they do not have 
terms that depend on primary variables Tq ). 

Comparing (I165l) - (ll66p and (I168p - (I169I) . we can see that all terms with temporal deriva- 
tives of the parameter (up to the second order which requires tertiary constraints, as 
we argued in 0, El) and temporal derivatives of primary fields, in the part of the generator 
(I165l) - fll66p we considered, coincide with diffeomorphism transformation (I168p - (I169I) without 
any need for field dependent redefinition of the gauge parameters. 

If the transformations for some fields that were calculated using the same generator 
require a field dependent redefinition of gauge parameters to "correspond" to diffeomorphism 
transformations, then such a redefinition will obviously destroy transformations (11650 - 01660 . 
If that happens then either the Hamiltonian formulation is not correct (ordinary mistakes 
are possible in such long calculations or non-canonical change of variables were performed) 
or the method of restoration of the gauge invariance is incorrect or, perhaps, sensitive to a 
choice of constraints; in that case such a method cannot be called an algorithm. 

There are a few more contributions to the transformations of different fields that we can 
obtain using the part of the generator (11441) . even without knowledge of the coefficient- 
functions a" and fl" This part of (11441) . after substitution of expressions for secondary 
constraints fl23l) and fl24l) . is 



G' = -i» X ™ = -e° {-h™ + h 00 n 00 - h km n km ) - e k (h°° k + 2h 00 n 0k + 2h m0 n km ) . (170) 



It allows us to find some contributions to the transformations of h' 



Sh at3 = {G\ h aP ) 



(171) 



that for different components leads to: 



53 



Sh pq = (-h km 7c km ) - s k (2h m0 TT km ) , h pq ) = -e°h pq + e p h q0 + e q h p0 , (172) 

6h °p = {-i k (2/i Q0 7r ofc ) , h 0p } = e p h 00 , (173) 

Sh 00 = (h 00 n 00 ) , h 00 } = i°h 00 . (174) 

Using the known transformations for g a ^ under diffeomorphism invariance and the defi- 
nition of h a P ({TBI , one obtains 

SMffhT = h*e% + h vX e$ - (h^e x ) x (175) 
which for different components of h a ^ gives: 

5 diff h°° = h 00 e° + 2h 0p e° p - e°h™ - e p h™ - h 00 e p p , (176) 

S diff h 0p = h 00 e p + h 0m e p m + h pm e° m - e°h 0p - h 0p e m m - e m h% , (177) 



S diff h pq = h p0 e q + h q0 e p + h pm e q m + h qm e p m - h pq e° - h pq e m m - e°h pq - e m h p ^ . (178) 

Similar to transformations of Tq fll65p - fll66l) . all found in the (1172p - (1174p contributions, 
are present in (1176p - (1178p and cover all terms with temporal derivatives of the gauge pa- 
rameter So, field dependent redefinition of gauge parameters is also not needed. 

The same part of the generator (11701) can be used to find partial transformations of 7i a/3 

5ir a p = {G',7i af3 } , (179) 

which gives for components 

8-K pq = {e°h km 7T km , ir pq } = e°-K pq , (180) 

5n 0p = {e°h°i: - e k 2h m0 TT km , n 0p } = ~e% - e k n kp , (181) 
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fcroo = {-e° (h 00 7T 00 ) - e k (h™ + 2/i 00 7r ofe ) , 7r 00 } = -£°7r 00 + e k k - 2e k ir 0k . (182) 

Returning to our original redefinitions, (fl5l) and (|T8|) . we find that ix a p = £ ai g. Now we 
can compare (I180ft - (I182I) with the transformations of 7t a p under diffeomorphism by using 

SdiffKap = 5rfj//S a/ g (183) 

and redefinition (j4j), which allows us to express H a p in terms of T^p (transformation of T^g 
under diffeomorphism is known (11671) ). From (11831) . for the components of 7T a p, we obtain: 

- e m Tip q ^ m - (27rpo + T^) e° q - 7t pm e™ - (2ir q0 + F™ m ) e° p - 7T qm e™ , (184) 

r 1 m I yvm 

Odiff"0p — ^ >P0 n pm^fl T ^Op^m 

OdiffKQO — — /7r mO£,o ~ ^ ^00,0 — ^ ^00,fc + J- 00 £ ,m ~ ^00^,0 • {i-OO) 

Here, as in (I168ft -( fl~69l) . the solutions for T™ and Y*™ must also be substituted (but it 
is not needed at this stage). Compare (I180l) - (ll82ft with the diffeomorphism transformation 
(jl84p - (|186P ; we see that again all terms with temporal derivatives of the gauge parameter e M 
are exactly the same, including two contributions with mixed spatio-temporal derivatives. 

Note also that the last two transformations, ( 1176j) -( Tl~78l) and ( 1 180h - ( TT82T) . follow just 
from an explicit form of the secondary constraints and obviously cannot be affected (in the 
Castellani procedure) by a choice of tertiary constraints. The transformations (1 165ft - ( 1166]) 
were obtained using the explicit form of coefficient-functions for which a particular choice of 
tertiary constraints (1 153ft was used; and there is possibility that different choices can affect 
the result. Of course, from physical point of view, this is impossible and if such happens, 
it would be an indication of problems with a method of the restoration of gauge invariance. 
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So, let us consider the effect of a different choice of constraints; e.g. one which we discussed 
in Section 4. To be explicit, let us consider (JSTj) for which the Hamiltonian takes the form 



K = ^ + A\™ + A m X Z (187) 
For this choice, contribution (11551) is modified 

, n-nk 00 , -00-pO ^0-00 „k00. ( -\ qq\ 

+ 2r oo r fc + T o r oo - «o r o - a o T k ; (18°) 

but solutions for coefficient-functions are the same as those given by (I156p - (I157I) . The 
remaining contributions calculated earlier, for the original choice of constraints, are also 
the same. One can easily check another choice that we considered in Section 4, (I85p . and 
verify that it also preserves the parts of coefficient-functions proportional to Tq . These 
particular examples illustrate the general statement made in Section 4 that a choice of 
linear combinations of constraints (or their PBs algebra) should not affect physical results 
(should not, in particular, affect gauge invariance). 

The above examples of partial restoration of gauge invariance using the Castellani proce- 
dure, which were considered, lead to partial contributions to the transformations of all fields 
presented in the reduced Hamiltonian of affine-metric formulation of GR. All contributions 
found are exactly the same as the corresponding terms of diffeomorphism transformations; 
and there is no need for a field dependent redefinition of gauge parameters. Moreover, all 
these contributions are independent of a choice of tertiary constraints, as it should be, if we 
respect the concept that the gauge invariance is a unique characteristic of a gauge invariant 
theory. Of course, to make the final conclusion, full calculations must be performed and all 
contributions have to be found. These calculations are straightforward, but quite laborious; 
the results will be reported elsewhere. Our main goal here is to find the reason for the 
contradictory result of [1JJ about necessity to have "correspondence" with diffeomorphism 
invariance. Now, after we show that diffeomorphism invariance follows directly from the 



Castellani procedure; but when used in [ll] the HTZ ansatz does not give correct transfor- 
mations. One can draw the conclusion that the HTZ ansatz cannot be considered to be an 
algorithm for the restoration of the gauge invariance and the need of field dependent redef- 
initions found in ll] is just an artifact of the method used. In next Section we explicitly 



demonstrate the failure of HTZ and discuss the possible reasons for this failure. 
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VII. HENNEAUX-TEITELBOIM-ZANELLI ANSATZ 



In this Section we turn our attention to the HTZ approach [12[ . This approach was used 



m [11] where the gauge invariance obtained for the first order Einstein-Hilbert action differs 
from the results of the previous Section where we follow the Castellani procedure. According 
to ll|], the gauge invariance of first order formulation of GR is not a four-diffeomorphism 
and a field dependent redefinition of gauge parameters is needed to "correspond to diffeo- 
morphism invariance" . 



The appearance of such a "correspondence" in |ll| cannot be explained by a non-canonical 
change of variables, as it was in the case of ADM model 4|. We have the Hamiltonian of the 
same theory and one method, Castellani (see previous Section), leads to four-diffeomorphism 



and another, HTZ (see III]), does not. This is clearly a deficiency of HTZ approach and 
this is exactly what we want to investigate and discuss in this Section. 

There are two major differences between the HTZ approach and the Castellani algorithm. 
The first one is the use of a so-called extended formalism: instead of the total Hamiltonian (as 
in Castellani's) the extended formalism is the starting point of the HTZ approach, where all 
primary constraints are included in the Hamiltonian with Lagrange multipliers. The second 
one is so-called HTZ ansatz - a generator of the gauge transformation is assumed to be a 
linear combination of all first class constraints and all of them enter explicitly {l^ . contrary to 
Castellani's generator fll36p - (ll39l) where only primar y fi rst class constraints explicitly enter 
a generator. Moreover, the HTZ iterative procedure 12j to find field dependent coefficients 
in front of the constraints (see (I189p ). is started from the end of the constraint chains (e.g. 
start from tertiary constraints for the first order affine-metric GR). 

The question of a "total versus extended" Hamiltonian will not be discussed here. We 
only would like to mention that in the particular application of the HTZ approach the 
"gauge fixing" of the Lagrange multipliers of non-primary first class constraints is used. 
This effectively converts the extended Hamiltonian into a total Hamiltonian (see Eq. (4.1) 
of [3]). Moreover, different variations of the HTZ method were developed; and some of 



them are based entirely on the total Hamiltonian, as in Castellani procedure (e.g. [721 ]) 
without even mentioning the extended Hamiltonian; but they lead to the same results as 
HTZ approach. So, our interest is the HTZ ansatz, which is also the essential part of all 



modifications of this approach (e.g. see Eq. (5) of 72|). In 11[ the gauge transformations 
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were obtained "using a method very similar to the method of HTZ" . We will consider the 
HTZ ansatz following the original paper [3] (the reprint of this paper can also be found in 
Sections 3.2-3.3 of the book [7l|). 

According to the HTZ ansatz [3] for a system with first class constraints, the generator 
is simply a linear combination of all first class constraints. As in the Castellani algorithm, 
all second class constraints should be eliminated in the preliminary step and the Dirac 
brackets should be calculated jl2j ]; this is exactly what we did in Section 2. In the case of 
;he Hamiltonian with tertiary constraints, the HTZ generator is given by (see Eq. (4.2a) of 



G = a«^i + ^ 2 ^ 2 + • (189) 

We slightly adjust the HTZ notation to make it more transparent for a covariant theory. 
(p^ are the first class constraints of different generations i = 1,2,3 (primary, secondary, and 
tertiary) and a Mi are functions of the canonical variables and inexpressible velocities (T^n in 



the considered formulation) that are iteratively defined from equation (see Eq. (4.2b) of 12] 
for i > 2) 



where C Vi ^ and V v ^ are structure functions in the PBs of the primary constraints with 
the rest of constraints {C Ui/3 ^) and in PB of any constraints with the canonical Hamiltonian 



{^,M = J2 C ^^ > ( 191 ) 

i<s 



i<s+l 



is a short notation for (see Eq. (3.4c) of [l2j) 



Dt ~ dt + ^dT^ + TU °°dt 00 + "■ (193) 
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According to [12|], one has to start from the equation of the highest order and "without 
loss of generality" take a M3 = e M , i.e. functions (gauge parameters), which are independent 
of the canonical variables that leads for % — 3 to 

— + {e^ h c } + {e», aUj Kr - ^ ^ c J; = °- ( 194 ) 

i>2 j>i 

We would like to emphasize that the independence of the gauge parameters of the canon- 
ical variables is the starting point of this iterative procedure; and without this assumption 
one will face the problem of solving variational equations instead of algebraic ones. Despite 
the computational problems with field dependent parameters, there are no a priori criteria 



of what possible dependence should be assumed. Let us, as the authors of [12J suggested 



and as it was done in 



ll| (see Eq. (145) of [ll|), take a M3 = e^; i.e. as independent of the 



phase-space variables function of spacetime coordinates. 

This equation, (11941) . can be simplified as the result of independence of of the canonical 
variables: = J^; and both PBs in (11941) are zero. Because the HTZ ansatz (11891) and 
structure functions (119ip -( fT92l) explicitly depend on a choice of constraints, we have to 
specify our choice from outset (note that in previous Section, using the Castellani algorithm, 
we were able to perform some calculations without referring to the explicit form of the 
tertiary constraints). Partial restoration of gauge invariance using the HTZ approach for 
first order GR was discussed in [11[ for a particular choice of tertiary constraints. So, we also 
consider the same combinations, (I58I) and (1591) (the same choice was used in previous Section 
which led directly to four-diffeomorphism when the Castellani procedure was applied). For 
these constraints the corresponding canonical Hamiltonian is 

i 1,0m 

H c = -n o xl° + ^ + ^t° + A\™ + A- X Z (195) 

The PBs among the chosen constraints are the simplest, {x°°) r °°} = 0, and primary 
constraints have zero PBs with all secondary and tertiary constraints making all C Vx p = 
(see (11911) ). which leads to a simple form of equation (I194p 

^-a^V^-e»V u ^=0. (196) 

The explicit form of the structure functions V v ^ (I192p for this choice of tertiary constraints 
is 
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1 h° 
— V ms = — — - 

which allows to solve equation (11961) . which is algebraic with respect to a 1 " 1 



K Z = ~^ V" = -7^ V p ™* = 5™ V^» = 0, (197) 



a 



a o 2 = -h 00 ^- + h°°e v V v ^ , (198) 

- h 0p ^- + ^£% 3 ° 3 - £ V V V P \ (199) 

Using (1781 - ( IHTl) . the remaining structure functions V u p3 can be easily found (note that 
they have only dependence on canonical variables h 0a ), which makes a p2 to be functions 
of the gauge parameters and variables h 0a . This immediately allows us to find the partial 
transformations of h a P using part of the generator fl 1 8 9 [) with a pl and a M2 (already found) 



Sh 



a 3 



{h* G} = {h* ... + a>»tf + e% 00 } = JL; (... + a^xl° + ^%°°) ■ (200) 



This is especially simple for the component h 0a , as the corresponding momenta are present 
only in the secondary constraints, which lead to 



5y 00 

5/l oo = a ^^_ = a o^oO ; 

07T UU 



£ 00 

5h 0p = fl M2_^_ = aP2/l 00_ 



After substitution of explicit form of V^ 3 into (I198p - (I199I) . we obtain: 



(201) 
(202) 



,02 



~h 



00 



de° . 

h e ■ 

dt ' m 



h 0r ' 



h 00 



— £ 



,0m 



h 00 



h 00 



„■ - = s p - h 0p i° + h 00 e pm ( — ] e° - e pm e° m - 



,P2 



1 



(203) 
(204) 



e m + 



h 0r. 



h 00 ,m h 00 



- ^S'Z, + h 0p ( -t- 



+ h 0p 



h 0m 



h 0m h 0p 
h 00 £ ' m 



Equation (12031) is exactly the same as Eq. (146) of [llj; but (I204[) has different signs in 
front of a few terms, as compared to the similar Eq. (147) of 11] . These transformations, 
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(l20ip - (j2Q2p . are not transformations under four-diffeomorphism (see (11 761) - (11771) in previous 
Section) and at most the so-called "correspondence" can be found by a field dependent 
redefinition of the gauge parameters. Note that such a "redefinition" is not needed when 



the Castellani a 
was found in 



Hi: 



gorithm is employed. The "appropriate" redefinition of gauge parameters 



£ ° = (205) 



h 0i 

* = ? ~ ( 206 ) 



Substitution of these expressions into (1203}) gives 

6h oo = h oog> Q + 2h o P ^ _ £O h oo _ eh oo _ h oo^ (2Q7) 

As we mentioned above, our (I204p has a few different signs compared to the similar Eq. 
(147) of Q; this can probably explain why only a transformation Sh 00 was provided in [ljj . 
Substitution of (12051) and (12061) into (T20¥|) gives 



8h°* = h 00 ^ + h 0m £% + h^i% - eti* - h^C m - Ch% . (208) 



According to the author of [ll| (I203j) "is IDENTICAL [Capital letters are ours] with 
diffeomorphism invariance ... IF we substitute" f!205j) and (I206p . i.e. perform field dependent 
redefinition of gauge parameters, which by the initial assumption of the HTZ ansatz should 
be independent of the phase-space variables. 

What is the significance of such a "correspondence", especially if there is a different 
method that leads directly to four-diffeomorphism? Firstly, by the transformations produced 
by the HTZ ansatz are (12031) - (I204p and their derivation was based on the assumption (that 
was used in the course of derivation) that the gauge parameters are field independent. 
So, a field dependent redefinition of the gauge parameters is just a manipulation, not a 
derivation; and this manipulation is in contradiction with what was used to derive the 
gauge transformations. Secondly, why is this particular redefinition of parameters chosen? 
If the gauge invariance of some theory is not known a priori, then it is meaningless to seek 
a correspondence to this unknown gauge invariance precisely because it is unknown. In 
such a case what should be called the gauge invariance if such manipulations are allowed? 
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Gauge invariance is a unique and very important characteristic of a theory, and neither a 
Hamiltonian formulation nor methods of restoration that lead to such ambiguities (or the 
need for such manipulations) can be accepted. Returning to the discussion in Section 5, 
it is quite obvious that the HTZ method cannot give a covariant result for this choice of 
tertiary constraints because these constraints, (j58p and (1591) . have different dimensions, as 
do the corresponding gauge parameters. In the Castellani procedure the choice of linear 
combinations of tertiary constraints is irrelevant and any combination leads to the same 
unique gauge transformation, whereas in the HTZ method, at least some of the combinations 
of tertiary constraints definitely cannot lead to a covariant result. This gives a limitation 
on the application of this method and imposes severe restrictions on the possible operations 
with non-primary first class constraints. 

Let us take a different combination of tertiary constraints which gives H" of (11871) . one 
that we used in previous Section where we showed that in this case the Castellani procedure 
produced the same gauge transformations as using any other combination. Let us return 
to (I196p . which is the first step of the iterative procedure, and its form is the same for any 
choice of tertiary constraints (only structure functions will be different). Now we consider 
the following linear combinations of tertiary constraints 

1 1,0m 

f oo _ _L r oo , « oo ~oo _ oo / 2n q^ 

~~ ^00 faOO m ' 'm — 'm ■ l^ UJ J 

In terms of these constraints the canonical Hamiltonian is 

tt _ -pv ..,,00 _,_ ~00 _,_ /lO^OO , Am. O0 fOirA 

ti c — -J-ooX^ + r o + A Xo + A X m - [ 2W ) 

In this case the structure functions V U ^ A of the HTZ ansatz become (as {H c , x°®} = — t®°+ 
terms proportional to x°u) 



Kf = (211) 



and the solution to equation (1 196[) is 



« M2 = -— + e v V/ A \ (212) 



~dt 

Even without specifying the structure functions V v ^ , which are more complicated compared 
to the previous choice of tertiary constraints, some contributions to the transformations can 
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be calculated using the general expression (jSJ) 



5h* = {h*G} = {h*... + ^ 7 ff + ..\ 



at s**- (213) 

This partial contribution to the transformations leads exactly to the same result as in 
(1180p - fll82p (for this choice of constraints, (12091) . equation (12131) equals to (11791) in previous 
Section), which are part of four-diffeomorphism without any need for redefinition of gauge 
parameters, in full correspondence with the original assumption of their independence of 
phase-space variables. 

This is a clear demonstration of the sensitivity of the HTZ ansatz to the choice of tertiary 
constraints. It cannot be considered an algorithm if such ambiguities are possible. We have 
already discussed this in Section 4 (see (1901) ) using the arguments that invariance of the 
corresponding Lagrangian cannot depend on a choice of tertiary constraints. 

The simple examples, (I20ip - (12041) and fl 2 1 3 j) . are related to calculation of a small part of 
the full generator and transformations that it produces, more discrepancies might appear 
if full calculations are performed. If the HTZ ansatz for each choice of linear combinations 
of tertiary first class constraints gives different transformations, it cannot be considered 
a reliable method of finding gauge invariance, which must be unique characteristics of a 
system. The only possibility to reconcile the ambiguities of the HTZ ansatz with a unique 
gauge invariance of a system is an existence of one particular choice of tertiary constraints 
and this choice must be specified in the HTZ approach. How can such a combination of 
tertiary constraints be found? It makes no sense to even try. Especially, since it is not clear 
at all whether it is necessary to do this if there exists the Castellani algorithm that allows 
us to work with any combination of tertiary constraints without affecting the unique gauge 
invariance of a system; and, at the same time, we can pick a combination of constraints with 
which it is easier to perform calculations. 

The possible existence of one special combination of constraints gives too strong restric- 
tion on the HTZ approach and is opposite to the advantages of the method that were stated 
in 12]. In particular, according to the authors 3], if the restriction that structure func- 
tions become structure constants is imposed (see Eq. (4.5b) [j^j]) then the generators can 
be written in the form given previously by Castellani. In reality, the opposite is true: only 
with exactly such a restriction, (12111) . the HTZ ansatz gives some meaningful results, con- 
trary to the Castellani algorithm, which is independent on a choice of linear combinations of 
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tertiary constraints, as it should be in a reliable algorithm. For completeness of discussion 
of the HTZ ansatz we provide our arguments on what restrictions might be imposed on a 
choice of constraints in the HTZ ansatz to recover the same invariance with the Castellani 
algorithm (which is free from such restrictions). These arguments are based on calculations 
that were performed in |4( where to restore the four-diffeomorphism in Dirac's Hamiltonian 
formulation of second order GR we used both methods: Castellani's and one based on the 
HTZ ansatz developed in 



721 ] . We obtained the same result, four-diffeomorphism with no 
need for any field dependent redefinition of parameters, by using both methods; and this, 
as we understood now, was an "accidental" result. The reason for this is a very special 
choice of secondary constraints, as in (11141) . that we used in |4j; and they were defined as 
the result of the PBs of primary constraints with the Hamiltonian. This is the only case 
where the HTZ ansatz might work. One realistic field-theoretical example considered by the 
authors of fl is Yang-Mil* theory, which a. 8 o .eaxis to known gangs transformations; it 
was also considered by Castellani [16J to illustrate the use of his algorithm. But in case of 
Yang-Mills theory, the secondary constraints were also defined as PBs of the corresponding 
primary constraints with the Hamiltonian. It seems to us that only with such a restriction 
on what we should call a secondary constraint the HTZ ansatz reproduces the same results 
as the Castellani algorithm. In the first order affine-metric formulation of GR, tertiary con- 
straints also appear. It is not clear how to find the "right" combination of constraints. The 
HTZ method might work if we define a "tertiary constraint" as everything that is produced 
as a time development of the corresponding secondary first class constraint. In contrast, 
according to the Dirac procedure at every step, when we consider a time development of 
known constraints, we have, first of all, to single out constraints which are already found, 
then the remaining part can be called a new constraint. In the formulation considered here, 
first order metric-affine GR, PBs of secondary first class constraints with the Hamiltonian 
give the expressions in which we can isolate terms proportional to a linear combination of 
secondary constraints. But this linear combination is not unique, and we do not have any 
prescription to find which linear combination is preferable. In addition, if all linear combina- 
tions of secondary constraints are equally good, why do we have to use a particular and very 
special linear combination of tertiary constraints? From our point of view, it is preferable 
to use the method which is free from these problems. 

Another important observation, based on our calculation in the second order EH action 



64 



using the Castellani or [72j methods is that the amount of calculation is the same in both 
cases; so there is no advantage in using the HTZ ansatz even if we can figure out what 
restrictions should be imposed for the general case with long chains of constraints. In 
contrast, our partial calculations in the previous Section show that the independence of 
the final result on which linear combinations of tertiary constraints are used, gives us an 
advantage as some of them can drastically simplify the calculations. 

The sensitivity of the HTZ ansatz and the problems related to it provide an additional 
illustration of importance of primary first class constraints in the Hamiltonian formulation 
that we discussed in Section 5. They are the true Masters of the Hamiltonian formulation 
and because the Castellani approach uses explicitly only primary constraints, it gives correct 
results and takes care of the possible redefinition of all non-primary constraints. But methods 
based on the HTZ ansatz, which starts from non-primary constraints, is an attempt to 
interchange the roles of constraints and this is the reason for its failure. The HTZ ansatz 
has some similarities with the non-canonical change of variables in the ADM formulation: 
both treat primary constraints as unimportant and both emphasize the role of non-primary 
constraints, either by starting the iterative procedure from them or by completely fixing one 
particular combination of secondary constraints that are reflected in giving the special names: 
"Hamiltonian" and "spatial diffeomorphism" constraints. Such names are in contradiction 
with the covariance of General Relativity and with the Dirac conjecture that all first class 
constraints are responsible for gauge invariance. For example, among eight constraints of 
the second order formulation of GR in four-dimensional spacetime, if we separate three of 
them and call them a "spatial diffeomorphism" constraint, for what invariance are the five 
remaining constraints responsible? 



The origin of the failure of the HTZ ansatz is made 



especia. 



ly clear from another article 



[73j |. which was published almost simultaneously with [12| . In [73] the Castellani procedure 
for construction of a generator of gauge transformations is considered, and the authors 
remark that "the problem is complicated by the fact that the chain algorithm [the name 
used for the Castellani algorithm f ll36H139p ] strongly depends on the representation of the 
primary first-class constraints surface that is adopted. More precisely, if the primary first- 
class constraints are 



'a l 



0, a 1 = l,...,m 1 (214) 
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and can each be taken as the head of a chain in [our (11361) ]. it is in general not true that the 
equivalent constraints 



4>' ai =M b a l(q,p) ( f )bl =0 (215) 

also lead to consistent solution of [our (1136111391) ]. So, algorithm is not invariant under [our 
(12151) ]. An algorithm ... that is invariant under redefinition [our (I215p ] has been recently 
proposed in [our jl2 |]." 

So, according to [73j the HTZ ansatz is built on the assumption of invariance under (I215P 
while the Castellani algorithm is not. 

Of course, the Hamiltonian formulation "strongly depends" on the primary first class 
constraints; but this is not a "problem". It is just the nature of the Hamiltonian formulation 
that is reflected in the Castellani algorithm 14 . As we discussed in Section 5, primary first 
class constrains are special due to their dual nature: they are constraints and, at the same 
time, they are phase-space variables. As constraints, they allow redefinition (12 1 5[) without 
any restriction. If primary constraints are pure momenta (as it happens in most cases), 
then (I215p would be also the redefinition of phase-space variables and it must be canonical! 
So this is not a surprise that an "algorithm" which allows an arbitrary redefinition (12151) . 
leads to loss of a unique gauge invariance because in general such a redefinition does not 
correspond to a canonical transformation; and consequently, the new total Hamiltonian will 
not be equivalent to the original one, as well as its original gauge invariance will not be 
restored without some field dependent redefinition of gauge parameters (at best). So any 
"algorithm" that is invariant under arbitrary (I215p will unavoidably lead to the loss of gauge 
invariance for majority of choices of M%} {q,p)- 

Calculations with the Hamiltonian of GR are involved and to illustrate the devastat- 



ing effect of (I215p the simple Hamiltonian formu 
Electodynamics for which the total Hamiltonian {5 



ation can be considered, e.g. Maxwell 

a is 



H T = i°7r - A°d k n k - ^n k n k + i (d k A m - d m A k ) (d k A m - d m A k ) 



14 This kind of problem is similar to what is stated in the first line of "Golden Oldie" 70|: "The general 
coordinate invariance underlying the theory of relativity creates basic problems in the analysis of the 
dynamics of the gravitational field" . ADM variables is the solution of the problem with covariance. 
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One can take, for example, the following redefinition 



7Tq = A iY k n 

and try to consider the time development of this new primary constraint, closure of the 
Dirac procedure, algebra of constraints and the restoration of gauge invariance in this case. 

VIII. CONCLUSION 

In this paper we presented a detailed derivation of the Hamiltonian for the first order 
affine-metric formulation of GR, including restoration of its unique gauge invariance, four- 
dimensional diffeomorphism, using the Castellani algorithm, and we demonstrated that four- 
diffeomorphism can be lost as the result of a non-canonical change of variables or by using 
methods of restoration which are sensitive to a choice of linear combinations of non-primary 
first class constraints. These results are based on mathematical derivations and as such do 
not need any interpretation or discussion. They can be disproved by indication of mistake(s) 
or must be accepted. One additional, "conventional", option is just to ignore them saying 
that it is the well-known fact that "the canonical treatment breaks the symmetry between 
space and time in general relativity and the resulting algebra of constraints is not the algebra 



of four diffeomorphism" 42j and, because of this, only by some unjustified manipulations 
the "correspondence to diffeomorphism invariance" can be accomplished. 

We address this article to the readers who make their judgement based on the results, 
not on the correspondence of results to "conventional wisdom". Such readers, as well as 
"conventional" ones, do not need a long conclusion to make their minds and, because of 
this, we stop our discussion here. 
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APPENDIX A: SOLVING THE SECOND CLASS CONSTRAINTS 

We outline the last step of the Hamiltonian reduction, the elimination of two additional 
pairs of the phase-space variables (n° m , E k m ) and (n™,r™). The part of the canonical 
Hamiltonian, H c , with terms proportional to and r? (two last lines of (|22p ). the only 
source of contributions into secondary second class constraints fl27|) and (1281) . is 

h c (£L, rL) = -h 00 ^ k 0m + 2h*° (r^ p £™ - i^e™) - h km 2^ 0m - 2/#e*. (ai) 
+h km (2r^rL - iy*, - r p kq ri P )-2h k \ 00 TT m +2h km (n 0p ri m - 2n 0k ri q )-h™r k pq +2h pk r 

The first line of (lAlj) can be presented in the following form 
where 

D m — h m + h q it qm + h q T qm — h T p mp . (A3) 

Note that is not a traceless combination. 

The secondary constraint fl2"T]) . can be obtained from the variation of (lAlj) (using 
the fundamental PB {j: k m ,U°P} from ([TOD) 

X TT Sm 
frOm _ fj-tOm tt \ _ U£lc _ v 0m _ tOOy™ , n 0m k r>0n / A /i \ 

^fc - j^fc > ^c) - - J^jf- ~Xk ~ h L 0fc + U k ~ D - x U n • l A4 J 

This constraint obviously has a non-zero PB with the primary constraint and this 
pair of constraints (n° m , Xom) * s °f second class and the corresponding pair of variables, 
(n° m , Sp m ) , can be eliminated. Solving = from (1A4I) for E^. and substituting the 
pair 
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n? = o, S£ = - ^-^c) (A5) 

into (|22il gives the next reduction with the following change in the canonical Hamiltonian 



H c , 



H c (XL) = H c (SL from Eq.QM)) = ~f^Df + ^^^r. (A6) 
Separating in Z)° a contributions proportional to T k m we write 

A°n = ^ + (rj m ) (A7) 

where 

£>™ = h% + h^ qm , (A8) 

and 

An (r£J = w°r£n - h ko r? mp . (A9) 

Note that the trace of Dj^ (r^J is zero and for (IA6I) ; we obtain 

jT \ j- n0« nQ» i * ^ nQ" nQa ^ nQa nQfr frfc ~\ _jL nOa (r k ^ nO b fr fc ~\ 

^c(X 0m ) =- — D b D a +JJZT[hOO IJn D * ~h00 Vb Da V-irnJ-^Ub qm) U a [} gm) • 

(A10) 

The first two terms in (1A101) correspond to the second line of ( 13T1) and the last two terms 
of (lAlOj) must be combined with the second line of (lAlj) . Now the part of the canonical 



Hamiltonian that depends on T k can be written as 



qui 



V (rk \ _ km ( r>-pg yp _ pP p? _ "np "ng \ _ j^qkym I A i i \ 

^c^V/ ~ V fcm 1 kp 1 mq 1 kq 1 mp) ly m 1 qk- K 1 ^ 11 ) 

In the quadratic part of ( 1A11I) , the combination e km naturally arises and , in the part 
linear in T^,, can be written in manifestly symmetric form 

hqo hko 

Dt = ~2h k H Qm + h» - w ,D» k - — (A12) 
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+ (V 7r 00 + 2/i**7Tqp - h% + j^c*) 5« m + (V°7r 00 + 2hPH Qp - h% + ^D° c ^j 5 h ri 



where is defined as (IA8I) (or f[53j) in the main text). This part of if c , flAHI) . leads to 
the secondary second class constraints 



iff = {Hf, H c } = = X f = 0. (A13) 

and the last pair of phase-space variables can be eliminated: IT^ = and the solution of 
Performing the variation 



flAT3|) for Performing the variation in flAT3l) we obtain 



2e , T « 9 _ e ky TL _ e kz Y y kx + 51 {e km T\ m - e^) + 51 {e km r km - e kz T%) = Df . (A14) 
The way to solve (IA14D for Tf is similar to the Einstein proof of equivalence of the first 

h n 

and second order formulations of metric GR [jj (see also [5(, Appendix A). The solution of 



flA14j) is based on the subsequent elimination of "traces". Contracting (1A14I) with 5* we 



obtain 

km-pz _ kz-pP _ \ jjxz (A~\^\ 
e 1 km e L kp — y) _ I x ' K 1 ^ 10 ) 

The left-hand side of (1A15I) is exactly the combination which appears in brackets with a 
Kronecker delta in equation ( 1A14I) . This allows us to write ( 1A14I) as 

2e y*Ti q - e k yri x - e kz Y\ x = Df - ^jD;y - ^-^f = Df. (A16) 

This combination, D^ z , is also used in the reduced Hamiltonian (13T|) ; and its explicit form 
is written in ( |34l) . 

We still have to eliminate the trace T xq in (1A16I) . Contracting ( 1A16I) with h yz we find 



1 

2 (D~- 2) 

and equation flA16j) takes the final form 



rs* = om oA ^' ( A17 ) 
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The solution of (1A18I) can be found using the Einstein permutation, as it was done in 



and 



(Appendix A). We need to have three free indices in the same position. To obtain 



this, we contract (1A18K with e nx to obtain 



- e nx e ky T z kx - e nx e kz T y kx = e nx D y x z . (A19) 

Now we perform a permutation in indices n, y and z and add the resulting equations in 
the following order (nyz) + (zny) — (yzn). That gives us 

e nx e kz T y kx = -- (e nx D y x z + e zx D™ v - e yx D zn ^j . (A20) 
Contracting (IA20j) with h an hb z we obtain the solution for Y y ba1 

r L = ~\ (h z D yz + h an Dl y - h an h bz e yx b zn ^) . (A21) 

The last pair of phase-space variables (Yl yz ,T x z ^ can be eliminated (the last Hamiltonian 
reduction). 

Substitution of ( 1A21I) into the part of the canonical Hamiltonian ( 1A11I) gives the third 



line of (EH) that is expressed in terms of the combination D yz . Note that in (1A11I) we have 
D yz and by using our redefinition ( 1A16I) can be also written in terms of D yz 



Df = Df - D»>81 - Dl z b y x . (A22) 
This completes the Hamiltonian reduction of the afhne-metric formulation of GR. 
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